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1.  INTRODUCTION 


Breast  cancer  is  incurable  upon  metastasis  to  distant  organs,  and  metastasis  to  axillary  lymph  nodes  is 
regarded  as  a  critical  prognostic  factor  for  future  recurrence  and  survival.  Understanding  the 
epidemiology  and  biology  of  metastasis  could  lead  to  better  stratification  of  recurrence  risk.  We 
proposed  to  study  genes  related  to  epithelial-mesenchymal  transition  (EMT),  invoking  the  hypotheses 
that  EMT  may  explain  the  ability  of  tumor  cells  to  form  metastatic  lesions  and  that  these  genes  are 
regulated  via  DNA  methylation.  It  is  hypothesized  that  tumor  cells  co-opt  the  EMT  program  to 
transiently  acquire  properties  generally  reserved  for  mesenchymal  cells,  namely  the  ability  to  detach  and 
migrate.  The  objectives  of  this  project  were  to  interrogate  the  protein  expression  and  promoter 
methylation  of  six  EMT-related  genes:  E-cadherin,  N-cadherin,  Vimentin,  Twistl,  RelB,  and  SATBl. 
Protein  expression  has  been  measured  via  immunohistochemistry  (IHC)  in  breast  tumor  tissue  and 
promoter  methylation  will  be  measured  using  DNA  derived  from  these  tumor  samples.  Protein 
expression  and  methylation  status  will  be  correlated  with  lymph  node  metastasis  at  diagnosis,  time  to 
metastatic  recurrence,  and  disease-free  survival.  Effect  modification  by  tumor  grade,  hormone  receptor 
status,  and  HER2  status  will  also  be  investigated. 

2.  KEYWORDS 

Breast  cancer;  Metastasis;  Epithelial-Mesenchymal  Transition;  Prognosis;  Molecular  Epidemiology; 
Methylation;  Immunohistochemistry 


3.  OVERAUU  PROJECT  SUMMARY 
Training  Plan 

Tasks  1/5  -  All  predoctoral  program  requirements  have  been  completed  and  dissertation  defense  is 
scheduled  for  April  2015. 

Task  2  -  Trainee  has  regularly  attended  and  participated  in  journal  clubs  (Cancer  Prevention  and 
Epigenetics  research  groups),  work  in  progress  meetings,  and  other  relevant  Institute  seminars.  Trainee 
attended  several  national  conferences,  including  the  American  Association  for  Cancer  Research  (AACR) 
Annual  Meeting  and  the  San  Antonio  Breast  Cancer  Symposium. 

Task  3  -  Trainee  has  conducted  several  additional  molecular  epidemiology  research  projects  focusing  on 
molecular  and  genetic  factors  relating  to  lymph  node  metastasis,  recurrence,  and  survival.  In  2012, 
trainee  participated  in  a  week-long  summer  course  in  survival  analysis  (Survival  Analysis  Applied  to 
Epidemiologic  and  Medical  Data,  University  of  Michigan  School  of  Public  Health  Graduate  Summer 
Session  in  Epidemiology). 
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Task  4  -  Trainee  has  worked  with  a  pathologist  at  RPCI  and  a  breast  biology  researcher  at  the 
University  at  Buffalo  (Dr.  Patricia  Masso-Welch)  to  identify  and  implement  an  appropriate  plan  for 
interpreting  immunohistochemical  stains  and  is  currently  completing  this  work. 


Summary  of  Results,  Progress,  and  Accomplishments: 

A.  An  analysis  of  the  relationship  between  tumor  size  and  lymph  node  metastasis  by  tumor 
subtype  in  breast  cancer  patients  of  African  and  European  ancestry  enrolled  in  the  Women’s 
Circle  of  Health  Study  (WCHS)  was  presented  at  the  2012  AACR  Annual  Meeting.  We 
found  that  European-American  (EA)  women  with  small  tumors  (<2  cm)  were  at  decreased 
likelihood  of  being  lymph  node  positive  at  diagnosis.  A  similar  trend  was  observed  for 
African-American  (AA)  women,  though  the  association  was  not  statistically  significant. 
Eurther,  we  found  that  ER  negativity  was  associated  with  decreased  risk  of  node  positivity 
(OR=0.41,  95%  Cl  0.20-0.84)  among  AA  women  with  large  tumors,  but  with 
nonsignificantly  increased  risk  among  AA  women  with  small  tumors  (OR=1.90,  95%  Cl 
0.92-3.91).  Partly  because  pathology  data  was  unavailable  for  all  participants,  this  analysis 
was  limited  by  sample  size,  particularly  for  subgroup  stratification,  and  may  have  therefore 
been  underpowered  to  detect  significant  associations.  Because  recruitment  of  AA  participants 
in  the  WCHS  has  been  ongoing  and  additional  participants  now  have  tumor  size  and  receptor 
status,  we  plan  to  revisit  this  project  in  the  near  future  [Appendix  1]. 

B.  In  a  second  project,  we  examined  the  relationships  between  seven  single  nucleotide 
polymorphisms  (SNPs)  in  the  BRMSl  and  SIPAl  genes  and  lymph  node  status,  tumor 
characteristics,  overall  survival,  and  recurrence-free  survival  in  a  cohort  of  859  women 
diagnosed  with  invasive  breast  cancer,  who  were  enrolled  in  the  Data  Bank  and 
BioRepository  at  Roswell  Park  Cancer  Institute.  These  results  were  presented  at  the  2011 
AACR  Annual  Meeting  and  at  the  2011  San  Antonio  Breast  Cancer  Symposium  [Appendices 
2  and  3]. 

The  manuscript  was  published  in  Breast  Cancer  Research  and  Treatment  in  2013  [Appendix 
4].  We  found  that  lymph  node  positive  tumors  were  less  likely  among  patients  with  the 
SIPAl  rs3741378  variant  genotype,  and  more  likely  among  patients  heterozygous  for  the 
BRMSl  rs  1052566  variant  (Table  2).  Having  the  variant  genotype  of  SIPAl  rs7 894763  was 
associated  with  an  increased  risk  of  high  grade  tumors  (Table  3).  Table  4  shows  associations 
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between  the  SNPs  and  tumor  subtype.  The  variant  genotype  of  BRMSl  rs31 16068  was 
associated  with  an  increased  risk  of  having  the  luminal  B  or  HER2-enriched  tumor  subtypes, 
while  the  BRMSl  rs  1052566  variant  was  associated  with  a  reduced  risk  of  the  luminal  B 
tumor  subtype.  The  variant  genotypes  of  SIPAl  rs746429  and  rs2306364  were  associated 
with  reduced  risk  of  the  triple  negative  subtype.  We  did  not  observe  any  significant 
associations  with  survival  or  recurrence  (Table  5).  Finally,  to  assess  the  effects  of  these  SNPs 
together,  we  created  a  summary  risk  allele  score  (Table  6).  We  found  that  having  8  or  more 
risk  alleles  was  associated  with  significantly  increased  risk  of  lymph  node  positive  tumor, 
and  that  overall,  there  was  a  dose-response  relationship  between  the  number  of  risk  alleles 
and  likelihood  of  node  positivity  (Ptrend  =  0.002).  There  were  no  significant  associations 
between  the  summary  score  and  tumor  grade  or  the  survival  outcomes,  however. 

C.  A  third  project  investigates  the  effect  of  polymorphisms  in  12  metastasis-related  genes  on  the 
risk  of  breast  cancer,  stratified  by  lymph  node  status  and  estrogen  receptor  status,  in  A  A  and 
EA  women  enrolled  in  the  WCHS.  Using  the  adaptive  rank  truncated  product  method  of 
pathway  analysis,  we  found  that  variants  in  the  CDHl  and  SIPAl  genes  were  significantly 
associated  with  risk  of  lymph  node  positive  and  ER  negative  breast  cancer,  respectively,  in 
AA  women.  In  EA  women,  we  identified  SNPs  in  the  MTA2,  SATBl,  KISSl,  SNAIl, 

CD82,  NMEl,  and  CTNNBl  genes  as  being  potentially  important  markers  of  lymph  node  or 
estrogen  receptor  status.  These  results  were  submitted  as  an  abstract  to  the  2015  AACR 
Annual  Meeting  [Appendix  5]. 

D.  A  fourth  project  proposes  to  analyze  26  genetic  variants  in  several  EMT-related  genes 
(YAPl,  AREG,  CDH2,  FOXMl,  SNAIl,  and  REEB)  in  relation  to  breast  tumor 
characteristics,  lymph  node  status  at  diagnosis,  and  recurrence  in  the  Pathways  Study,  a  large 
cohort  study  of  breast  cancer  survivorship.  Genotyping  has  been  completed  and  data  analysis 
is  ongoing. 

Current  Objectives: 

A.  Finalize  manuscript  describing  the  results  presented  in  Appendix  5  and  submit  for 
publication.  This  project  forms  the  second  chapter  of  the  trainee’s  dissertation. 
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B.  Continue  work  on  ongoing  projects  (analysis  of  tumor  size  and  nodal  status  relationship  in 
WCHS  and  of  genetic  variants  in  the  Pathways  Study).  Plan  to  submit  these  papers  for 
publication  in  the  coming  year. 


Research  Plan 

Task  1  -  Interpretation  of  IHC  assays  is  ongoing,  using  a  combination  of  manual  scoring  and  the 
positive  pixel  count  algorithm  provided  by  Aperio.  Analysis  of  the  IHC  data  will  be  the  trainee’s  third 
dissertation  chapter. 

Task  2  -  We  initially  received  DNA  from  458  participants  with  tumor  tissue  from  the  Pathology  core 
facility  at  Roswell  Park.  We  examined  these  samples  for  quality  and  quantity  by  several  different 
methods,  and  found  that  the  samples  contained  insufficient  DNA  for  methylation  analysis.  We  therefore 
requested  new  FFPE  cores  for  DNA  preparation.  Primary  tumor  and  matched  metastatic  lymph  node 
cores  for  a  subset  of  the  patient  population  have  been  received.  We  are  currently  exploring  options  to 
obtain  funding  for  methylation  analysis.  Because  these  matched  samples  are  a  unique  resource,  we  plan 
to  measure  methylation  of  the  proposed  six  genes,  as  well  as  investigate  methylation  of  additional  loci. 
Tasks  3/4  -  Analysis  of  the  IHC  data  will  first  be  conducted  on  the  matched  primary  tumor  and 
metastatic  lymph  node  samples.  This  analysis  will  be  part  of  the  trainee’s  dissertation.  Following 
successful  defense,  the  remainder  of  the  primary  tumor  data  will  be  analyzed  for  relationships  with 
recurrence  and  survival.  We  expect  to  present  these  data  at  upcoming  conferences  and  to  have  the  paper 
published  by  early  2016. 

Task  5  -  Over  the  grant  period,  the  trainee  has  presented  several  posters  at  major  conferences  and  has 
presented  data  at  seminars  and  work-in-progress  meetings  at  Roswell  Park  and  the  University  at  Buffalo. 
The  trainee  has  also  actively  participated  in  journal  clubs  to  aid  in  the  development  of  good  scientific 
communication  and  presentation  skills. 

Summary  of  Results,  Progress,  and  Accomplishments: 

A.  All  tumor  tissue  was  received  as  expected  and  immunohistochemistry  for  the  6  proteins  was 
completed.  Stained  tissue  microarray  slides  were  scanned  using  the  Aperio  ScanScope  in  the 
Pathology  core  facility.  The  interpretation  of  the  immunohistochemistry  panel  is  nearing 
completion. 
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Current  Objectives: 

A.  Complete  interpretation  of  immunohistochemistry  and  analyze  data  with  respect  to  tumor 
characteristics  and  prognosis,  with  planned  publication  by  early  next  year. 

B.  Investigate  and  apply  for  funding  to  examine  methylation  status  at  these  6  loci  as  well  as 
other  metastasis-related  loci. 

4.  KEY  RESEARCH  ACCOMPLISHMENTS 

Nothing  to  report. 

5.  CONCLUSION 

Results  from  the  genetic  association  studies  conducted  by  the  trainee  indicate  that  inherited  variants  in 
metastasis-related  genes  may  affect  tumor  characteristics,  in  particular  the  propensity  to  form  metastases 
to  the  axillary  lymph  nodes.  Because  presence  of  lymph  node  metastases  is  a  critical  prognostic  factor 
for  breast  cancer  patients,  it  is  crucial  to  understand  the  mechanisms  by  which  tumor  cells  gain 
metastatic  potential.  Our  results  need  to  be  replicated  in  additional,  larger  patient  populations;  we  may 
be  able  to  leverage  existing,  publicly  available  genotyping  datasets  to  further  investigate  the  role  of 
inherited  genetic  variation  on  metastatic  potential  as  well  as  to  replicate  our  findings.  We  anticipate 
completing  analysis  of  the  immunohistochemistry  data  by  the  end  of  the  year.  Here,  too,  we  may  be  able 
to  utilize  existing  gene  expression  repositories  to  replicate  our  findings.  Finally,  we  are  presently 
exploring  options  to  acquire  funding  for  more  extensive  methylation  analyses  than  we  previously 
proposed,  using  analytes  from  the  primary  tumor  and  matched  metastatic  cores  we  received.  In  addition, 
we  can  leverage  this  resource  to  explore  other  genomic  variation,  such  as  copy  number  variants.  Our 
hope  is  that  this  information  can  lead  to  improved  prognostic  stratification  for  breast  cancer  patients. 

6.  PUBLICATIONS,  ABSTRACTS,  AND  PRESENTATIONS 

a.  Manuscripts 

1)  Lay  press:  Nothing  to  report 

2)  Peer-reviewed  Scientific  Journals: 

i.  Roberts  MR,  Hong  CC,  Edge  SB,  Yao  S,  Bshara  W,  Higgins  MI,  Freudenheim  JL, 
and  Ambrosone  CB.  Case-only  analyses  of  the  associations  between  polymorphisms 
in  the  metastasis-modifying  genes  BRMSl  and  SIPAl  and  breast  tumor 
characteristics,  lymph  node  metastasis,  and  survival.  Breast  Cancer  Research  and 
Treatment;  2013  Jun;  139:3:873-85.  DOI:  10.1007/sl0549-013-2601-3.  PMID: 
23771732. 

3)  Invited  Articles:  Nothing  to  report 
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4)  Abstracts: 

i.  Roberts  M,  Hong  CC,  Edge  S,  Yao  S,  Nesline  M,  and  Ambrosone  CB. 
Polymorphisms  in  metastasis  suppressor  genes  (BRMSl  and  SIPAl):  Breast  tumor 
characteristics  and  lymph  node  metastasis  [abstract].  In:  Proceedings  of  the  102nd 
Annual  Meeting  of  the  American  Association  for  Cancer  Research;  201 1  Apr  2-6; 
Orlando,  FL.  Philadelphia  (PA):  AACR;  Cancer  Res  2011;71(8  Suppl): Abstract  nr 
5595.  DOI:  10. 1 158/1538-7445. AM201 1-5595. 

ii.  Roberts  MR,  Hong  CC,  Edge  SB,  Yao  S,  Nesline  M,  and  Ambrosone  CB.  Single 
nucleotide  polymorphisms  in  the  BRMSl  and  SIPAl  metastasis  suppressor  genes  as 
prognostic  markers  in  breast  cancer  patients.  In:  Proceedings  of  the  34^'’  Annual 
CTRC-AACR  San  Antonio  Breast  Cancer  Symposium;  2011  Dec  6-10;  San  Antonio, 
TX.  Philadelphia  (PA):  AACR;  Cancer  Res  2011;71(24  Suppl): Abstract  nr  Pl-09-06. 
DOI:  I0.II58/0008-5472.SABCSII-PI-09-06. 

iii.  Roberts  M,  Bandera  EV,  Hwang  H,  Ciupak  G,  Zirpoli  GR,  Yao  S,  Pawlish  K,  Davis 
W,  Jandorf  E,  Bovbjerg  DH,  and  Ambrosone  CB.  Tumor  size  and  lymph  node 
metastases  in  African-American  and  European- American  women  with  breast  cancer 
[abstract].  In:  Proceedings  of  the  I03rd  Annual  Meeting  of  the  American  Association 
for  Cancer  Research;  2012  Mar  31-Apr  4;  Chicago,  IE.  Philadelphia  (PA):  AACR; 
Cancer  Res  2012;72(8  Suppl):Abstract  nr  3593.  DOI:1538-7445.AM2012-3593. 

iv.  Roberts  MR,  Sucheston-Campbell  EE,  Zirpoli  GR,  Bandera  EV,  Ambrosone  CB,  and 
Yao  S.  Single  nucleotide  variants  in  metastasis-related  genes  are  associated  with 
breast  cancer  risk,  by  lymph  node  involvement  and  ER  status,  in  women  with 
European  and  African  ancestry  [abstract].  To  be  presented  at  the  106*  Annual 
Meeting  of  the  American  Association  for  Cancer  Research;  2015  Apr  18-22. 

b.  Presentations  made  during  the  past  year:  Nothing  to  report. 

7.  INVENTIONS,  PATENTS  AND  LICENSES 

Nothing  to  report. 


8.  REPORTABLE  OUTCOMES 

•  One  manuscript  published  in  Breast  Cancer  Research  and  Treatment  [Appendix  4] 

•  Two  abstracts  presented  at  AACR  Annual  Meetings  [Appendices  1  and  3] 

•  One  abstract  presented  at  the  San  Antonio  Breast  Cancer  Symposium  [Appendix  2] 

•  One  abstract  to  be  presented  in  April  2015  [Appendix  5] 

9.  OTHER  ACHIEVEMENTS 

•  Trainee  will  complete  her  doctoral  degree  in  May  2015 

•  Trainee  was  awarded  an  annual  research  award  of  $1500  (the  Saxon  Graham  Research  Award) 
by  the  University  at  Buffalo  Department  of  Social  and  Preventive  Medicine  in  May  2011 
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3.  Roberts  M,  Bandera  EV,  Hwang  H,  Ciupak  G,  Zirpoli  GR,  Yao  S,  Pawlish  K,  Davis  W,  Jandorf  E, 
Bovbjerg  DH,  and  Ambrosone  CB.  Tumor  size  and  lymph  node  metastases  in  African-American  and 
European- American  women  with  breast  cancer  [abstract].  In:  Proceedings  of  the  I03rd  Annual 
Meeting  of  the  American  Association  for  Cancer  Research;  2012  Mar  3I-Apr  4;  Chicago,  IE. 
Philadelphia  (PA):  AACR;  Cancer  Res  2012;72(8  Suppl): Abstract  nr  3593.  DOI:1538- 
7445.AM2012-3593. 

4.  Roberts  MR,  Hong  CC,  Edge  SB,  Yao  S,  Bshara  W,  Higgins  Ml,  Freudenheim  JE,  and  Ambrosone 
CB.  Case-only  analyses  of  the  associations  between  polymorphisms  in  the  metastasis-modifying 
genes  BRMSl  and  SIPAl  and  breast  tumor  characteristics,  lymph  node  metastasis,  and  survival. 
Breast  Cancer  Research  and  Treatment;  2013  Jun;  139:3:873-85.  DOI:  I0.I007/sI0549-0I3-260I-3. 
PMID:  23771732. 

11.  TRAINING  OR  FELLOWSHIP  AWARDS 

Training  activities:  Trainee  participated  in  and  presented  at  weekly  journal  clubs  and  work-in-progress 
meetings.  In  learning  new  statistical  approaches  for  the  analysis  of  genetic  data,  the  trainee  worked  with 
Lara  Sucheston-Campbell,  an  epidemiologist  in  the  Department  of  Cancer  Prevention  at  Roswell  Park 
Cancer  Institute.  Trainee  completed  a  one- week  course  in  survival  analysis  to  further  develop  data 
analysis  skills.  To  learn  how  to  interpret  immunohistochemistry  results  and  develop  a  scoring  system, 
the  trainee  worked  with  Patricia  Masso-Welch,  an  Associate  Professor  in  the  Department  of 
Biotechnical  and  Clinical  Laboratory  Sciences  at  the  University  at  Buffalo,  who  specializes  in  mammary 
gland  biology. 

Professional  activities:  Trainee  attended  and  presented  research  at  several  international  conferences. 
Trainee  also  continues  to  attend  numerous  regular  seminars,  including  weekly  seminar  series  at  Roswell 
Park  and  the  University  at  Buffalo,  and  monthly  meetings  specific  to  breast  cancer  research  at  Roswell 
Park. 
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12.  APPENDICES 


Appendix  1 :  Abstract  and  poster  entitled  “Tumor  size  and  lymph  node  metastases  in  African- 
American  and  European-American  women  with  hreast  cancer” 


Michelle  R.  Roberts^  Elisa  V.  Bandera^,  Helena  Hwang\  Gregory  Ciupak\  Gary  Zirpoli\  Song  Yao\ 
Karen  Pawlish^,  Warren  Davis',  Lina  Jandorf"',  Dana  H.  Bovbjerg^,  Christine  B.  Ambrosone' 

'Roswell  Park  Cancer  Institute,  Buffalo,  NY 

^University  of  Medicine  and  Dentistry  of  New  Jersey,  New  Brunswick,  NJ 
^New  Jersey  Department  of  Health  and  Senior  Services,  Trenton,  NJ 
"'Mount  Sinai  School  of  Medicine,  New  York,  NY 
^University  of  Pittsburgh  Cancer  Institute,  Pittsburgh,  PA 

Few  studies  have  evaluated  lymph  node  metastasis  in  African-American  (AA)  women  with  breast 
cancer,  who  are  more  likely  to  be  diagnosed  at  an  advanced  stage  and  with  lymph  node  positive  tumors. 
Likelihood  of  nodal  involvement  increases  with  tumor  size,  although  recent  data  have  indicated  that  this 
may  not  be  true  for  AA  breast  cancer  patients  and  patients  with  basal-like  tumors.  Nodal  metastases  are 
also  more  likely  in  premenopausal  AA  patients  compared  to  either  premenopausal  European-American 
(LA)  patients  or  postmenopausal  AA  and  EA  patients.  We  examined  risk  factors  for  lymph  node 
metastasis  at  breast  cancer  diagnosis  in  AA  and  EA  women,  and  investigated  the  contributions  of  race, 
tumor  subtype,  and  menopausal  status  to  the  tumor  size-lymph  node  metastasis  relationship.  This 
analysis  included  805  women  diagnosed  with  primary,  incident  breast  cancer  enrolled  in  the  Women's 
Circle  of  Health  Study,  a  case-control  study  of  AA  and  EA  breast  cancer  patients  and  healthy  women. 
Cases  were  identified  using  hospital-based  ascertainment  in  New  York  City  hospitals  with  high  referral 
patterns  for  AA  women  and  through  population-based  ascertainment  in  New  Jersey  using  the  State 
Cancer  Registry.  Eligible  cases  were  self-identified  AA  and  EA  women  age  20-75  with  no  previous 
history  of  cancer  other  than  nonmelanoma  skin  cancer.  In-person  interviews  were  conducted  and 
consent  to  review  pathology  reports  was  obtained.  Tumor  size  was  categorized  as  tumors  2  cm  or  less 
(small  tumors)  and  tumors  greater  than  2  cm  (large  tumors).  Unconditional  logistic  regression  was  used 
to  estimate  odds  ratios  (OR)  and  95%  confidence  intervals  (Cl).  AAs  with  small  tumors  were  more 
likely  to  be  node  positive  compared  to  EAs  (OR=1.24,  95%  Cl  0.81-1.88)  while  among  patients  with 
large  tumors,  AAs  were  less  likely  to  be  node  positive  (OR=0.83,  95%  Cl  0.49-1.41).  When  grouped  by 
race  and  tumor  subtype,  we  found  that  the  triple  negative  subtype  was  associated  with  a  decreased  risk 
of  nodal  metastases  among  EA  women  with  small  tumors  (OR=0.17,  95%  Cl  0.04-0.79)  and  a 
nonsignificantly  decreased  risk  among  AA  women  with  large  tumors,  using  the  luminal  A  subtype  as  the 
referent  group.  Associations  were  null  in  EA  women  with  large  tumors  and  AA  women  with  small 
tumors.  When  grouped  by  race  and  ER  status,  ER  negativity  was  associated  with  a  decreased  risk  of 
nodal  metastases  among  AA  women  with  large  tumors  (OR=0.41,  95%  Cl  0.20-0.84),  while  AA  women 
with  small  tumors  were  at  increased  risk  (OR=1.90,  95%  Cl  0.92-3.91).  Our  data  suggest  an  effect  of 
race  and  tumor  subtype  on  the  relationship  between  tumor  size  and  likelihood  of  lymph  node  metastases. 
Tumor  size  appears  to  affect  lymph  node  metastasis  differently  by  race,  a  mechanism  that  is  modified  by 
tumor  biology.  Our  findings  support  the  hypothesis  that  in  AA  breast  cancer  patients,  large  tumors  may 
not  be  more  likely  to  give  rise  to  metastatic  lymph  nodes. 

See  attached  poster  at  end  of  report. 
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Appendix  2:  Abstract  and  poster  entitled  “Polymorphisms  in  metastasis  suppressor  genes  (BRMSl 
and  SIPAl):  Breast  tumor  characteristics  and  lymph  node  metastasis” 


Michelle  R.  Roberts^  Chi-Chen  Hong\  Stephen  B.  Edge^,  Song  Yao\  Mary  Nesline\  Christine  B. 
Ambrosone^ 

'Department  of  Cancer  Prevention  and  Control,  Roswell  Park  Cancer  Institute,  Buffalo,  NY 
^Department  of  Surgical  Oncology,  Roswell  Park  Cancer  Institute,  Buffalo,  NY 

Introduction:  BRMSl  and  SIPAl  function  as  metastasis  suppressors,  but  few  studies  have  examined 
metastasis  suppressor  gene  polymorphisms  in  breast  cancer.  Axillary  lymph  node  (LN)  metastases  and 
tumor  characteristics  predict  aggressiveness  but  correlate  imperfectly  with  likelihood  of  metastatic 
relapse.  BRMSl  regulates  transcription  through  NF-kB  pathways;  protein  expression  has  been 
correlated  with  ER/PR  negative,  HER2  positive  tumors.  SIPAl  affects  extracellular  matrix  gene 
expression,  and  polymorphisms  have  been  associated  with  EN  metastases  and  ER/PR  negative  tumors. 
Identifying  polymorphisms  that  affect  metastasis  may  help  to  better  recognize  patients  who  require 
aggressive  adjuvant  therapy.  We  assessed  associations  between  SNPs  in  BRMSl  and  SIPAl  and  EN 
metastases,  tumor  grade,  and  ER/PR/HER2  status  in  breast  cancer  patients. 

Methods:  We  included  1,015  newly  diagnosed  breast  cancer  patients  who  received  surgery  at  Roswell 
Park  Cancer  Institute  (RPCI)  and  participated  in  the  DataBank  and  BioRepository  shared  facility. 
Participants  completed  an  epidemiologic  questionnaire  and  provided  a  blood  sample  prior  to  surgery  or 
other  treatment.  Clinical  and  pathologic  data  were  linked  to  de-identified  participant  data  in  the  DBBR 
database.  BRMSl  (rsl  1537993  and  rs31 16068)  and  SIPAl  (rs75894763)  SNPs  were  genotyped  through 
RPCFs  Genomics  shared  facility  using  Sequenom®  iPLEX  Gold  assays.  Logistic  regression  was  used 
to  estimate  odds  ratios  and  95%  confidence  intervals. 

Results:  Tumors  were  more  frequently  node  positive  among  never-users  of  hormone  replacement 
therapy,  and  node  positive  tumors  were  more  likely  to  be  high  grade,  ER/PR  negative,  and  HER2 
positive.  Node  positive  disease  was  less  likely  among  patients  heterozygous  for  BRMSl  rs31 16068 
(OR=0.76,  95%  CI=0.55-1.05)  and  SIPAl  rs75894763  (OR=0.42,  95%  CI=0.16-1.09).  BRMSl 
rsl  1537993  was  not  associated  with  lymph  node  metastases.  The  SIPAl  rs75894763  variant  allele  was 
also  associated  with  reduced  risk  of  ER  and  PR  negative  tumors.  HER2  positive  tumors  were  more 
likely  among  patients  homozygous  for  BRMSl  rs31 16068  (OR=2.38,  95%  CI=1.15-4.94).  No 
significant  associations  with  tumor  grade  were  observed. 

Conclusions:  Preliminary  data  indicate  that  polymorphisms  in  BRMSl  and  SIPAl  are  associated  with 
ER,  PR,  and  HER2  tumor  status,  and  with  reduced  risk  of  LN  involvement  among  breast  cancer  patients. 
Future  research  to  evaluate  these  and  other  genetic  variants  in  metastasis  suppressor  genes  in  relation  to 
recurrence  and  survival  is  necessary  to  better  understand  the  biology  of  metastasis. 

See  attached  poster  at  end  of  report. 
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Appendix  3  :  Abstract  and  poster  entitled  “Single  nucleotide  polymorphisms  in  the  BRMSl  and 
SIPAl  metastasis  suppressor  genes  as  prognostic  markers  in  breast  cancer  patients” 


Michelle  R.  Roberts^  Chi-Chen  Hong\  Stephen  B.  Edge^,  Song  Yao^  Mary  Nesline^  Christine  B. 
Ambrosone^ 

'Department  of  Cancer  Prevention  and  Control,  Roswell  Park  Cancer  Institute,  Buffalo,  NY 
^Department  of  Surgical  Oncology,  Roswell  Park  Cancer  Institute,  Buffalo,  NY 

Introduction:  Single  nucleotide  polymorphisms  (SNPs)  in  the  metastasis  suppressors  BRMSl  and  SIPAl 
may  affect  metastatic  efficiency.  BRMSl  affects  apoptosis,  colonization,  cell  adhesion,  and  invasive 
potential.  Loss  of  BRMSl  expression  has  been  correlated  with  younger  age  at  diagnosis  and  reduced 
survival  time  in  patients  with  progesterone -receptor  (PR)  negative,  HER2-positive  tumors.  SIPAl  affects 
extracellular  matrix  gene  expression  and  cell  adhesion,  and  while  SNPs  have  been  associated  with  node 
positive,  estrogen-receptor  (ER)/PR  negative  tumors,  evidence  for  a  relationship  with  survival  has  been 
conflicting.  Identifying  SNPs  that  affect  risk  of  recurrence  and  survival  may  improve  the  ascertainment 
of  patients  who  require  aggressive  adjuvant  therapy  following  a  diagnosis  of  breast  cancer.  We 
evaluated  associations  between  seven  BRMSl  and  SIPAl  SNPs  and  recurrence  and  survival  in  patients 
with  primary  breast  cancer. 

Methods:  We  identified  1,015  incident  breast  cancer  patients  who  received  surgery  at  Roswell  Park 
Cancer  Institute  (RPCI)  and  participated  in  the  DataBank  and  BioRepository  (DBBR)  resource. 
Participants  completed  an  epidemiologic  questionnaire  and  provided  a  blood  sample  prior  to  surgery  or 
other  treatment.  Clinical  and  pathologic  data  were  linked  to  de-identified  participant  data  in  the  DBBR 
database.  SNPs  in  BRMSl  (rsl  1537993,  rs31 16068,  and  rsl052566)  and  SIPAl  (rs75894763,  rs746429, 
rs3741378,  and  rs2306364)  were  genotyped  by  RPCI’s  Genomics  facility  using  Sequenom®  iPLEX 
Gold  and  Taqman®  real-time  PCR  assays.  Cox  proportional  hazards  regression  was  used  to  estimate 
hazard  ratios  and  95%  confidence  intervals. 

Results:  The  median  follow-up  time  was  33  months,  and  49  deaths  and  42  recurrences  occurred. 

Tumors  were  more  likely  to  be  larger,  node  positive,  ER/PR  negative,  and  high  grade  among  patients 
who  experienced  a  recurrence  or  death.  Recurrence  was  less  likely  in  older  patients  and  those  with 
higher  body  mass  index,  although  the  latter  association  was  nonsignificant  (p=0.06).  Patients  with  at 
least  one  variant  allele  of  the  BRMSl  rs31 16068  genotype  experienced  shorter  overall  survival 
compared  to  patients  with  the  homozygous  common  genotype  (HR=2.05,  95%  Cl  1.15-3.63,  rs3 116068 
AG-i-AA  compared  to  GG).  The  remaining  SNPs  were  not  associated  with  overall  survival,  and  none  of 
the  SNPs  were  associated  with  recurrence. 

Conclusions:  In  our  data,  the  variant  allele  of  rs31 16068  was  more  common  among  women  whose  breast 
cancer  was  node  negative  and  HER2 -positive,  compared  to  those  with  the  common  rs31 16068  allele. 

The  rs31 16068  variant  allele  is  also  associated  with  poorer  survival.  While  our  findings  do  not  support  a 
role  for  common  SNPs  in  the  SIPAl  gene  in  breast  cancer  prognosis,  BRMSl  rs31 16068  may  be  a  useful 
prognostic  biomarker.  Euture  goals  are  to  examine  additional  SNPs  in  BRMSl  and  other  metastasis- 
related  genes  in  a  larger,  racially  diverse  population. 

See  attached  poster  at  end  of  report. 
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Appendix  4\  Manuscript  entitled  “Case-only  analyses  of  the  associations  between  polymorphisms  in 
the  metastasis-modifying  genes  BRMSl  and  SIPAl  and  breast  tumor  characteristics,  lymph  node 
metastasis,  and  survival” 

See  attached  manuscript  at  end  of  report. 
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Appendix  5:  Abstract  entitled  “Single  nucleotide  variants  in  metastasis-related  genes  are  associated 
with  breast  cancer  risk,  by  lymph  node  involvement  and  ER  status,  in  women  with  European  and 
African  ancestry” 


Michelle  R.  Roberts \  Lara  E.  Sucheston-Campbell\  Gary  R.  Zirpoli\  Elisa  V.  Bandera^,  Christine  B. 
Ambrosone\  Song  Yao^ 

'Department  of  Cancer  Prevention  and  Control,  Roswell  Park  Cancer  Institute,  Buffalo,  NY 
^Rutgers  Cancer  Institute  of  New  Jersey,  New  Brunswick,  NJ 

Inherited  genetic  variation  may  partially  explain  inter-patient  variability  in  prognosis  by  influencing 
lymph  node  involvement  and  estrogen  receptor  (ER)  status  in  breast  cancer  patients,  which  may  differ 
by  ancestral  background.  We  examined  154  tagging  single  nucleotide  polymorphisms  (SNPs)  in  12 
metastasis-related  genes  (BRMSl,  CDHl,  CD82,  CTNNBl,  KISSl,  MTAl,  MTA2,  MTA3,  NMEl, 
SATBl,  SIPAl,  SNAIl)  for  associations  with  risk  of  breast  cancer,  stratified  by  lymph  node  and  ER 
status. 

Genotyping  was  performed  in  2,67 1  European- American  (EA)  and  African-American  (AA)  women 
enrolled  in  the  Women’s  Circle  of  Health  case-control  study  (WCHS)  using  Illumina  GoldenCate 
assays.  Single-SNP  and  haplotype  associations  were  analyzed  using  logistic  regression.  Pathway 
analyses  were  conducted  using  the  adaptive  rank  truncated  product  (ARTP)  method,  with  p<0.10  as 
significant.  To  estimate  risk,  multiallelic  scores  were  created  using  the  SNPs  in  the  significant  gene(s). 
All  models  were  adjusted  for  age  and  ancestry;  multiallelic  score  models  also  included  demographic 
covariates.  P-values  were  corrected  using  the  false  discovery  rate  (EDR)  method. 

Single-SNP  and  haplotype  associations  were  not  significant  after  EDR  adjustment  at  p<0.05.  In  AA 
women,  significant  ARTP  gene-level  associations  included  CDHl  with  risk  of  lymph  node  positive 
breast  cancer  (p=0.10)  and  SIPAl  with  ER  negative  breast  cancer  in  both  case-control  (p=0.10)  and 
case-case  (p=0.09)  analyses.  Multiallelic  scores  computed  from  SNPs  in  CDHl  and  SIPAl  were 
associated  with  node  positive  (OR=1.13,  95%  Cl  1.07-1.19,  pfdr=0.0003)  and  ER  negative  (OR=1.16, 
95%  Cl  1.02-1.31,  pfdr=0.03)  breast  cancer,  respectively. 

In  EA  women,  MTA2  was  associated  with  overall  risk  of  breast  cancer  at  the  ARTP  gene-level 
(p=0.004),  regardless  of  ER  status,  and  with  node  negative  breast  cancer  (p=0.01).  SATBl  and  KISSl 
were  also  significant  in  ER  negative  (ARTP  gene-level  p=0.03)  and  node  negative  (ARTP  gene-level 
p=0.10)  analyses,  respectively.  Among  EA  lymph  node  positive  cases,  significant  ARTP  gene-level 
associations  were  observed  for  SNAIl  (p=0.10),  CD82  (p=0.05),  NMEl  (p=0.10),  and  CTNNBl 
(p=0.09).  The  SNAI1-CD82-NME1-CTNNB1  multiallelic  risk  score  was  associated  with  node  positive 
(OR=1.09,  95%  Cl  1.04-1.14,  pfdr=0.001)  and  the  MTA2-K1SS1  score  with  node  negative  breast  cancer 
(OR=1.18,  95%  Cl  1.08-1.29,  pfdr=0.002).  Stratified  by  ER  status,  the  MTA2  score  was  associated  with 
ER  positive  (OR=1.16,  95%  Cl  1.05-1.30,  pfdr=0.01)  and  the  MTA2-SATB1  score  with  ER  negative 
breast  cancer  (OR=1.12,  95%  Cl  1.05-1.20,  pfdr=0.003). 

Our  findings  suggest  that  genetic  variants  in  several  metastasis  genes  may  affect  risk  of  breast  cancer  by 
lymph  node  or  ER  status.  These  results  require  verification  in  larger  studies,  particularly  those  that  can 
evaluate  long-term  prognosis. 
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Background  and  Methods 

'  Lymph  node  status  is  an  important  predictor  of  prognosis,  and  likelihood  of  nodal  involvement 
increases  with  tumor  size. 

'  Previous  data  indicates  that  larger  tumors  may  not  predict  risk  of  nodal  metastases  among  African- 
American  patients  or  those  with  basal-like  tumors. 

'  The  purpose  of  this  analysis  was  to  investigate  the  contributions  of  race,  menopausal  status,  and 
tumor  subtype  to  the  tumor  size-lymph  node  metastasis  relationship. 

'  Women  with  breast  cancer  were  enrolled  in  the  Women's  Circle  of  Health  Study  (WCHS),  a  case- 
control  study  recruiting  hospital-based  cases  from  New  York  City  hospitals  and  population-based 
cases  through  the  New  Jersey  State  Cancer  Registry. 

'  Self-identified  African-American  and  European-American  women  age  20-75  with  no  previous  history 
of  cancer  other  than  nonmelanoma  skin  cancer  were  eligible  for  enrollment. 

'  This  analysis  was  limited  to  721  cases  diagnosed  with  stage  I,  2,  and  3  invasive  breast  cancer,  who 
were  enrolled  between  2003  and  201 1  and  have  available  lymph  node  status  and  tumor  size. 

'  Unconditional  logistic  regression  was  used  to  estimate  age-adjusted  odds  ratios  (OR)  and  95% 
confidence  intervals  (Cl)  for  the  likelihood  of  having  a  node  positive  tumor  at  diagnosis. 


Table  2.  Race  and  Tumor  Size  are  Associated  with  Increased  Likelihood 
of  Lymph  Node  Metastases 


Characteristic 

Node 

Positive 

(N=286) 

Node 

Negative 

(N=435) 

OR  (95%  Cl) 

Race 

European-American 

African-American 

124  (43.3) 
162  (56.6) 

209  (48.0) 
226  (52.0) 

1.00 

1.31  (0.97-1.78) 

Menopausal 

Status 

Postmenopausal 

Premenopausal 

121  (42.3) 
165  (57.7) 

249  (57.2) 
185  (42.5) 

1.00 

0.97  (0.62-1.51) 

Tumor 

Subtype 

Luminal  A 

Luminal  B 

HER2-enriched 
Triple  negative 

128  (58.2) 

32  (14.5) 

16  (7.3) 

44  (20.0) 

217  (62.2) 
35  (10.0) 
25  (7.2) 

72  (20.6) 

1.00 

1.29  (0.75-2.22) 
0.99  (0.50-1.96) 
0.88  (0.56-1.38) 

Tumor  Size 

<2cm 

>2cm 

137  (47.9) 
149  (52.1) 

324  (74.5) 
111  (25.5) 

1.00 

3.07  (2.22-4.24) 

Table  3.  Association  Between  Tumor  Size  and  Lymph  Node  Status  by  Race, 
Menopausal  Status,  and  Tumor  Subtype 


Likelihood  of  Node  Positive  Tumor  by  Race,  Menopausal  Status,  and  Tumor  Subtype, 
OR  (95%  Cl) 


Tumor  Size 

African-American 

European-American 

Premenopausal 

Postmenopausal 

<2cm 

1.00 

1.00 

1.00 

1.00 

>2cm 

2.38  (1.55-3.64) 

4.27(2.55-7.14) 

2.28  (1.46-3.56) 

4.21  (2.63-6.75) 

Tumor  Size 

Luminal  A 

Luminal  B 

HER2-Enriched 

Triple  Negative 

<2cm 

1.00 

1.00 

1.00 

1.00 

>2cm 

4.41  (2.67-7.29) 

2.1  1  (0.76-5.86) 

1.00  (0.23-4.25) 

2.48(1.14-5.43) 

Table  1.  Characteristics  of  WCHS  Participants 


African-American  (N 

=371) 

European-American  (N= 

=350) 

Characteristic 

Node  Positive 

Node  Negative 

Node  Positive 

Node  Negative 

(N=162) 

(N= 

209) 

(N 

=  124) 

(N=: 

226) 

Age  at  diagnosis,  mean  (SD) 

48.3 

(9.8) 

52.5 

(10.4)t 

50.1 

(10.0) 

54.3 

(lO.O)t 

BMI,  mean  (SD) 

31.3 

(6.7) 

31.4 

(6.9) 

27.0 

(5.9) 

26.8 

(5.2) 

Age  at  menarche,  mean  (SD) 

12.6 

(2.0) 

12.5 

(1.8) 

12.4 

(1.6) 

12.4 

(1.4) 

. 

Nulliparous 

33 

(20.4) 

30 

(14.4) 

is" 

(30.7) 

66 

(29.2) 

an  y 

Parous 

129 

(79.6) 

179 

(85.7) 

86 

(69.4) 

160 

(70.8) 

Menopausal 

Premenopausal 

91 

(56.2) 

81 

(38.9)t 

74 

(59.7) 

104 

(46.0)t 

status 

Postmenopausal 

71 

(43.8) 

127 

(61.1) 

50 

(40.3) 

122 

(54.0) 

Yes 

18 

(11.1) 

34 

(16.3) 

32 

(25.8) 

51 

(22.6) 

Family  history 

No 

144 

(88.9) 

175 

(83.7) 

92 

(74.2) 

175 

(77.4) 

Positive 

90 

(64.8) 

125 

(65.1) 

90" 

(82.6) 

163 

(78.7) 

ER  status 

Negative 

49 

(35.3) 

67 

(34.9) 

19 

(17.4) 

44 

(21.3) 

PR  status 

Positive 

73 

(51.8) 

104 

(55.3) 

81 

(74.3) 

143 

(71.1) 

Negative 

68 

(48.2) 

84 

(44.7) 

28 

(25.7) 

58 

(28.9) 

Positive 

28 

(21.2) 

35 

(19.9) 

23 

(23.2) 

26 

(13.8) 

HER2  status 

Negative 

104 

(78.8) 

141 

(80.1) 

76 

(76.8) 

162 

(86.2) 

Low 

13 

(8.3) 

28 

(14.3) 

is" 

(15.9) 

66 

(30.1)t 

Tumor  grade 

Moderate 

59 

(37.8) 

75 

(38.3) 

54 

(47.8) 

87 

(39.9) 

High 

84 

(53.9) 

93 

(47.5) 

41 

(36.3) 

65 

(29.8) 

s2  cm 

68 

(42.0) 

133 

(63.6)t 

69" 

(55.7) 

191 

(84.5)t 

>2  cm 

94 

(58.0) 

76 

(36.4) 

55 

(44.4) 

35 

(15.5) 

1 

4 

(2.5) 

131 

(63.0)t 

11 

(8.9) 

187 

(82.7)t 

Stage  at 
diagnosis 

2 

84 

(51.9) 

77 

(37.0) 

77 

(62.1) 

38 

(16.8) 

3 

74 

(45.7) 

0 

(0.0) 

36 

(29.0) 

1 

(0.4) 

tp<0.001;  ^p<0.0001;  p-values  from  t-test.  Chi-squared  or  Fisher's  exact  tests,  as  appropriate. 

One  stage  4  patient  is  included  in  the  category  for  African-American  women  with  stage  3,  node  positive  tumors. 
BMIsbody  mass  index:  ER=estrogen  receptor  PR=progesterone  receptor;  AA=African-American;  EA=European-American. 


'  African-American  race  and  tumor  size  >2cm  were  associated  with  increased  likelihood  of  lymph  node  metastases.  High  grade  tumors  were 
more  likely  to  be  node  positive  in  European-American  women  only  (Tables!  and  2). 

'  Tumor  size  >2cm  was  significantly  associated  with  increased  likelihood  of  nodal  metastases  regardless  of  race  or  menopausal  status.  By 
subtype,  tumor  size  was  significantly  associated  with  nodal  metastases  in  the  luminal  A  and  triple  negative  subtypes.  There  was  no  effect  of 
tumor  size  in  the  HER2-enriched  subtype,  however  (Table  3). 


Figure  1.  The  tumor  size  and  iymph  node  status 
association  differs  by  race  and  ER  status. 


Figure  2.  Associations  between  race,  menopausai  status, 
tumor  subtype  and  nodai  status,  stratified  by  tumor  size. 


Figure  3.  Effect  of  ER-  and  HER2+  status  on  iymph  node 
status,  stratified  by  race  and  tumor  size. 


metastases  among  European-Americans,  but 
not  African-Americans. 


M.  Roberts  is  a  DOD  Predoctoral  Training  Fellowship  recipient  (W8IXWH-1 1 -1-0024).  WCHS  is  supported  by  USAMRMC  (DAMD-1 7-01-1 -0334),  NCI  (ROI  CAI 00598),  BCRF.  and  a  gift  From  the  Philip  L  Hubbell  family.  The 
NJ  State  Cancer  Registry  participates  in  the  Center  For  Disease  Control  and  Prevention's  National  Program  of  Cancer  Registries  and  is  a  NCI  SEER  Expansion  Registry,  supported  under  cooperative  agreement  DP07- 
703,  awarded  to  the  NJ  Department  of  Health  and  Senior  Services.  Collection  of  NJ  cancer  incidence  data  is  supported  by  SEER  under  contract  NOI -PC-95001 -20.  The  authors  have  no  conflict  of  interest  to  declare. 


differ  between  subgroups  defined  by  tumors  <2cm 
and  >2cm. 

In  HER2-enriched  tumors  >2cm,  there  was  a 
suggestion  of  decreased  risk  of  nodai  metastases. 


negative  tumors  >2cm  were  iess  iikely  to  be  node 
positive. 

in  European-American  women,  HER2  positive 
tumors  <2cm  ,  but  not  >2cm,were  more  iikeiy  to  be 
node  positive. 
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Abstract  #  5595 


ed  proliferation-associated  1  (SIPAI)  genes  supprc 
ence  and  survival,  but  few  data  are  av 


Background  and  Methods 

•Breast  cancer  metastasis  suppressor  1  (BRMS1)  and  s 
formation. 

•Axillary  lymph  node  metastases  are  an 
the  impact  of  genetic  variation  on  nodal 
•7  SNPs  In  BRMS1  and  SIPA1  were  evaluated  in  relation  to  lymph  node  metastasis  and  tumor  ch 
incident,  primary  breast  cancer  at  RPCI  from  2003-2010  and  who  are  enrolled  in  the  DataBank  a 
•1,015  patients  were  genotyped  with  Sequenom®  MassARRAY  IPLEX  Gold  assays  (rs1 1537993,  i 
BiosystemsTaqMan®  real-time  PCR  assays  (rsl  052566,  rs746429,  rs3741378,and  rs2306364). 


i  metastatic 
lable  to  evaluate 


Results  and  Conclusions 

•BRMS1  rsl  0  52566  was  associated  with  a 


I  at  diagnosis  (OR=1.42;  95%  Cl  1.02-1.96). 
observed  for  participants  with  the  variant  allele  of  S/PAI  rs3741378 


racteristics  in  women  diagnosed  with 
id  BloRepository  shared  facility. 
>3116068,  and  rs75894763)  and  Applied 


Table  1 .  Demographic  and  Tumor 
Node  Metastases  at  Diagnosis. 

Characteristi 

cs  by  Presence 

of  Lymph 

Characteristic 

’,n(%) 

Node  Positive 
(N  =  248) 

Node  Negativ 
(N  =  731) 

<50 

97  (39.1) 

226  (30.9) 

Age  at 

51-65 

90  (36.6) 

315(43.1) 

diagnosis 

>66 

61  (24.6) 

190  (26.0) 

White 

172  (90.5) 

566  (94.3) 

Race 

Nonwhite 

18(9.5) 

34  (5.7) 

High  school  or  less 

80  (42.1) 

192  (32.6) 

Education 

Some  college  or  graduate 

110(57.9) 

397  (67.4) 

0.02 

Menopausal 

Premenopausal 

72  (38.5) 

207  (34.6) 

0  34 

sUtus 

Postmenopausal 

115(61.5) 

391  (65.4) 

Never 

140  (75.7) 

398  (67.7) 

HRT  use 

Ever 

45  (24.3) 

190  (32.3) 

0.04 

Nulliparous 

28  (15.0) 

119(19.9} 

Parity 

Parous 

159  (85.0) 

480  (80.1) 

0.16 

Positive 

186  (75.0) 

537  (79.2) 

ER  status 

Negative 

62  (25.0) 

141  (20.8) 

0.18 

Positive 

155  (62.5) 

474  (69.9) 

PR  status 

Negative 

93  (37.5) 

204  (30.1) 

0.04 

Positive 

44(17.7) 

67  (11.2) 

HER2  status 

Negative 

204  (82.3) 

533  (88.8) 

0.01 

Well 

47  (19.2) 

181  (28.8) 

Tumorgrade 

Moderate 

108  (44.1) 

274  (43.6) 

0.004 

Poor 

90  (36.7) 

173  (27.6) 

Tis,  Tmi ,  T1A(<  5mm) 

14(5.6) 

210  (28.9) 

Tumor  size 

TIB,  T1C(>  5-20mm) 

116(46.8) 

417(57.3) 

<0.0001 

T2,  T3  (>  20mm) 

118(47.6) 

101  (13.9) 

HRT=hormone  replacement  therapy;  ER=estrogen  receptor;  PR=progesterone  receptor. 

'  Race,  education,  menopausal  status,  HRT  use,  and  parity  were  available  for  823  participants; 
age  at  diagnosis,  ER,  PR,  HER2  status,  tumor  grade,  and  tumor  size  were  available  for  1,015 
participants. 

^  p-values  from  Chi-squared  and  Fisher’s  exact  tests. 


lase  in  risk  of  lymph  node  me 

•Decreased,  borderline  statistically  significant  risk  of  lymph  node  positivity  wa 
{OR=0.70;95%  Cl  0.48-1.02). 

■The  variant  allele  of  SI  PA  1  rs3741378  has  previously  been  associated  with  estrogen  and  progesterone  positivity,  but  not  with  lymph  node  status. 

•BRMS1  rs3116068  was  associated  with  an  increased  risk  of  having  the  HER2  overexpressing  subtype  (OR=2.32;  95%  Cl  1.14-4. 73),  defined  as  a  tumor 
w hich  is  estrogen  and  progesterone  receptor  negative  and  HER2  positive. 

•These  findings  suggest  that  the  BRMS1  and  StPA  /  genes  may  affect  metastasis  to  regional  lymph  nodes  aswellastumor  receptor  status. 

•Future  studies  will  examine  genetic  variation  across  these  and  other  metastasis  suppressor  genes  to  confirm  these  findings. 


Table  2.  BRMS1  and  SIPA1  Single  Nucleotide  Polymorphism  Associations  with  Presence  of  Lymph  Node  Metastases,  Tumor  Receptor  Subtype,  and 
Tumor  Grade  at  Diagnosis. 


Genotype 

Odds  Ratio  for  Risk  of  Node 
Tumor  at  Diagnosis 

Positive 

Odds  Ratio  for  Risk  of  Luminal  B,  HER2  Overexpressing,  and 
Triple  Negative  Subtype  Versus  Luminal  A  Subtype  (n=604) 

Odds  Ratio  for  Risk  of  High  and  Moderate 
Grade  Versus  Low  Grade  Tumor  (n=234) 

positive 
n  (%) 

Node 
negative 
n  (%) 

Adjusted'  OR 
(95%  Cl) 

Adjusted*^  OR 
(95%  Cl) 

Adjusted^  OR  (95%  Cl) 
Luminal  B 
(n=76) 

Adjusted^ 
HER2  Ove 
(n 

OR  (95%  Cl) 
rexpressing 
=35) 

Adjusted^  OR  (95%  Cl) 
Triple  Negative 
{n=143) 

Adjusted''  OR  (95%  Cl) 
Moderate  Grade 
(n=386) 

Adjusted"  OR  (95%  Cl) 

H  igh  Grade 
(n=267) 

BRMS1 

AA 

128(51.8) 

366  (50.6) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

7811537993  (Leu67Leu) 

AG-fGG 

119(48.2) 

358  (49.4) 

0.92 

(0.67,  1.26) 

0.96 

(0.67,  1.44) 

1.27 

(0.78,  2.07) 

0.88 

(0.44,  1.78) 

0.78 

(0.53, 1.14) 

1.39 

(0.97, 1.98) 

1.17 

(0.75,  1.84) 

GG 

167  (68.2) 

444  (61.8) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

rs3116068  (3’  UTR) 

AG  +AA 

78  (31.9) 

275  (38.3) 

0.79 

(0.56,  1.11) 

0.69 

(0.46, 1.03) 

0.93 

(0.56,  1.55) 

2.32 

(1.14,  4.73) 

0.97 

(0.65,  1.45) 

1.05 

(0.73, 1.51) 

0.87 

(0.54,  1.40) 

GG 

104  (43.2) 

347  (49.5) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

rs1052566  (Ala273Val) 

AG  +AA 

137(56.8) 

354  (50.5) 

1.42 

(1.02, 1.96) 

1.50 

(1.02,  2.20) 

0.59 

(0.36,  0.96) 

0.60 

(0.29,  1.23) 

1.11 

(0.76,  1.64) 

0.90 

(0.63, 1.29) 

0.88 

(0.56,  1.40) 

SIPA1 

GG 

241  (98.0) 

691  (95.3) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

rs75894763  (Val621Val) 

AG 

5(2.0) 

34  (4.7) 

0.39 

(0.14,1.06) 

0.51 

(0.16, 1.62) 

1.43 

(0.53,  3.88) 

N/As 

0.29 

(0.07,  1.24) 

0.40 

(0.16,  0.96) 

1.13 

(0.43,  2.96) 

GG 

113(46.1) 

327  (46.2) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

rs746429  (Ala920Ala) 

GA-fAA 

132(53.8) 

381  (53.8) 

1.00 

(0.73, 1.38) 

1.10 

(0.75, 1.61) 

0.83 

(0.51,  1.36) 

0.54 

(0.26, 1.09) 

0.76 

(0.52,  1.11) 

0.73 

(0.51,  1.05) 

0.89 

(0.56, 1.41) 

CC 

190  (77.9) 

519(73.2) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

rs3741378  (Ser182Phe) 

TC  +  TT 

54  (22.1) 

190  (26.8) 

0.70 

(0.48, 1.02) 

0.72 

(0.46, 1.14) 

1.17 

(0.68,2.01) 

0.66 

(0.26, 1.64) 

1.20 

(0.79,  1.83) 

0.98 

(0.65, 1.48) 

1.36 

(0.82,  2.29) 

GG 

189  (77.5) 

526  (75.1) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

rs2306364  (Ala342Ala) 

AG  +AA 

55  (22.5) 

174(24.9) 

0.89 

(0.61, 1.31) 

0.84 

(0.54, 1.31) 

0.70 

(0.39, 1.28) 

0.77 

(0.33, 1.84) 

0.63 

(0.39, 1.02) 

1.29 

(0.85,  1.95) 

1.11 

(0.65, 1.91) 

Luminal  A:  ER  and/or  PR  positive,  HER2  negative;  Luminal  B:  ER  and/or  PR  positive,  HER2  positive;  HER2  overexpressing:  ER  and  PR  negative,  HER2  positive;  Triple  negative:  ER,  PR,  and  HER2  negative. 

Receptor  status  was  missing  for  162  participants. 

Unconditional  logistic  regression  was  used  to  estimate  odds  ratios  and  95%  confidence  intervals  for  risk  of  node  positive  tumors,  adjusting  for  ''age  at  diagnosis,  ER/PR  status,  tumor  size,  and  tumor  grade,  and  ^age  at 
diagnosis,  ER/PR  status,  tumor  size,  tumor  grade,  race,  education,  parity,  and  hormone  replacement  therapy  use. 

Multinomial  logistic  regression  (generalized  logit  model)  was  used  to  estimate  odds  ratios  and  95%  confidence  intervals  for  risk  of  luminal  B,  HER2  overexpressing,  or  triple  negative  tumor  subtype,  using  luminal  A  as  the 
comparison  group;  and  for  risk  of  high  grade  (poorly  differentiated  and  undifferentiated)  or  moderate  grade  tumors,  using  low  grade  (well-differentiated)  as  the  comparison  group. 

^Adjusted  models  included  age  at  diagnosis,  tumor  size,  and  lymph  node  status. 

^Adjusted  models  included  age  at  diagnosis,  tumor  size,  lymph  node  status,  and  ER/PR/HER2  status. 

^No  participants  had  both  heterozygous  genotype  and  HER2  overexpressing  subtype  tumors. 

OR=odds  ratio;  CNconfidence  interval. 

Supported  by  a  gift  from  the  Jayne  and  Phil  Hubbell  Family.  The  DBBR  and  Genomics  facility  are  RPCI  Cancer  Center  Support  Grant  shared 
resources  (P30  CA016056-32).  M.  Roberts  is  a  DOD  Predoctoral  Training  Fellowship  recipient  (BC100068)  and  C.  Ambrosone  is  funded  by  the 
Breast  Cancer  Research  Foundation.  The  authors  have  no  conflict  of  interest  to  declare. 
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Objectives  and  Methods 


't'Breast  cancer  metastasis  suppressor  1  (BRMS1)  and  signal-induced 
proliferation-associated  1  {SIPA1)  are  metastasis  suppressor  genes,  known  to 
suppress  metastatic  formation. 


*t*We  hypothesized  that  single  nucleotide  polymorphisms  (SNPs)  in  these  genes 
could  influence  recurrence  and  survival  in  women  with  breast  cancer. 


*:*Our  objective  was  to  evaluate  associations  between  7  SNPs  in  BRMS1  and 
SIPA1  and  time  to  recurrence  (TTR)  and  overall  survival  (OS). 


*MTR  —  time  from  diagnosis  to  date  of  first  recurrence  or  last  follow-up. 
*t*OS  —  time  from  diagnosis  to  date  of  death  or  last  follow-up. 


{'Participants  were  women  diagnosed  with  incident,  primary  breast  cancer  at 
RPCI  from  2003-2010  who  enrolled  in  the  DataBank  and  BioRepository  (DBBR) 
shared  facility. 


{'Genotyping  was  performed  using  Sequenom®  MassARRAY  iPLEX  Gold  assays 
(rs11537993,  rs3116068,  and  rs75894763)  and  Applied  BiosystemsTaqMan®  real¬ 
time  PCR  assays  (rs1052566,  rs746429,  rs3741378,  and  rs2306364). 


{'Pathology,  epidemiologic,  treatment,  and  follow-up  data  were  obtained  from  the 
DBBR  and  Tumor  Registry  databases  at  RPCI. 


't'Cox  proportional  hazards  regression  was  used  to  estimate  hazard  ratios  (HR) 
and  95%  confidence  intervals  (Cl)  for  relationships  between  the  SNPs,  TTR,  and 
OS  in  patients  with  stage  l-lll  breast  cancer  (N=850). 


't'Covariates  were  identified  using  log-rank  tests  of  significance. 


Results  and  Conclusions 


't'38  recurrences  and  44  deaths  occurred  during  follow-up  (median  follow-up 
time  was  33  months). 


't'There  were  no  statistically  significant  associations  with  time  to  recurrence. 


'Mn  general,  associations  with  overall  survival  were  null,  with  the  exception  of 
Si7Af5f  rs3116068  (HR=2.18,  95%  Cl  1.05-4.53,  for  heterozygous  compared  to 
homozygous  wildtype  genotype). 


't'SNPs  chosen  in  these  two  metastasis  suppressor  genes  do  not  appear  to 
correlate  with  recurrence  or  survival,  with  the  possible  exception  of  BRMS1 
rs31 16068. 


't'Our  analyses  were  limited  by  few  events  and  few  participants  with  the 
homozygous  variant  genotypes,  despite  a  relatively  large  sample  size. 


't'Longer  follow-up  times  are  likely  needed  to  further  investigate  SNPs  in  these 
genes  as  well  as  to  confirm  the  potential  association  of  BRMS1  rs31 16068  with 
overall  survival. 


Figure  1.  SNP  Associations  with  Time  to  Recurrence 
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Figure  2.  SNP  Associations  with  Overaii  Survivai 


Tabie  1.  Participant  Characteristics 


Participant  Characteristics,  n  (%) 

TTR  Events, 
n  (%) 

Censored, 

n  (%) 

OS  Events, 
n  (%) 

Censored, 
n  (%) 

<50 

19 

(50.0) 

257 

(31.7) 

10 

(22.7) 

266 

(33.0) 

Age  at  diagnosis 

51-65 

11 

(29.0) 

327 

(40.3) 

15 

(34.1) 

323 

(40.1) 

>66 

8 

(21.1) 

228 

(28.1) 

19 

(43.2) 

217 

(26.9) 

Race 

White 

28 

(93.3) 

606 

(93.4) 

31 

(88.6) 

603 

(93.6) 

Non-white 

2 

(6.7) 

43 

(6.6) 

4 

(11.4) 

41 

(6.4) 

Education 

High  school  or  less 

12 

(40.0) 

232 

(36.2) 

16 

(45.7) 

228 

(35.6) 

At  least  some  college 

18 

(60.0) 

409 

(63.8) 

19 

(54.3) 

408 

(64.2) 

<  25 

14 

(53.9)^ 

204 

(32.9) 

9 

(29.0) 

209 

(34.0) 

BMI 

25  -  29.9 

8 

(30.8) 

198 

(31.9) 

11 

(35.5) 

195 

(31.7) 

&  30 

4 

(15.4) 

218 

(35.2) 

11 

(35.5) 

211 

(34.3) 

Menopausal 

Premenopausal 

10 

(34.5) 

223 

(34.6) 

6 

(17.7)1 

227 

(35.5) 

status 

Postmenopausal 

19 

(65.5) 

422 

(65.4) 

28 

(82.4) 

413 

(64.5) 

HRT  use** 

Ever 

10 

(52.6) 

182 

(44.4) 

11 

(40.7) 

181 

(45.0) 

Never 

9 

(47.4) 

228 

(55.6) 

16 

(59.3) 

221 

(55.0) 

Parity 

Parous 

23 

(79.3) 

532 

(82.4) 

29 

(85.3) 

526 

(82.1) 

Nulliparous 

6 

(20.7) 

114 

(17.7) 

5 

(14.7) 

115 

(17.9) 

Family  history  of 

Yes 

4 

(13.3) 

139 

(21.4) 

6 

(17.1) 

137 

(21.3) 

breast  cancer 

No 

26 

(86.7) 

510 

(78.6) 

29 

(82.9) 

507 

(78.7) 

ER  status 

Positive 

20 

(54.1)2 

642 

(79.4) 

25 

(56.8)2 

637 

(79.4) 

Negative 

17 

(46.0) 

167 

20.6) 

19 

(43.2) 

165 

(20.6) 

PR  status 

Positive 

13 

(35.1)2 

559 

(69.1) 

14 

(31.8)2 

558 

(69.6) 

Negative 

24 

(64.9) 

250 

(30.9) 

30 

(68.2) 

244 

(30.4) 

HER2  status 

Positive 

4 

(10.8) 

105 

(13.1) 

8 

(18.2) 

101 

(12.8) 

Negative 

33 

(89.2) 

694 

(86.9) 

36 

(81.8) 

691 

(87.3) 

Lymph  node 

Positive 

16 

(42.1) 

230 

(28.4) 

16 

(42.9) 

228 

(28.3) 

status 

Negative 

22 

(57.9) 

580 

(71.6) 

24 

(57.1) 

578 

(71.7) 

Low 

(2.7)2 

215 

(27.0) 

4 

(9.3)2 

212 

(26.8) 

Tumor  grade 

Moderate 

15 

(40.5) 

356 

(44.7) 

15 

(34.9) 

356 

(45.0) 

High 

21 

(56.8) 

226 

(28.4) 

24 

(55.8) 

223 

(28.2) 

Tis,  Tmi,  T1 A  (<5mm) 

6 

(15.8)1 

95 

(11.7) 

3 

(6.8) 

98 

(12.2) 

Tumor  size 

TIB,  TIC  (&5-20mm) 

16 

(42.1) 

515 

(63.6) 

24 

(54.6) 

507 

(63.1) 

T2,  T3  {>20mm) 

16 

(42.1) 

200 

(24.7) 

17 

(38.6) 

199 

(24.8) 

Radiation 

Yes 

27 

(71.0) 

621 

(77.1) 

27 

(62.8)1 

621 

(77.6) 

No 

11 

(29.0) 

184 

(22.9) 

16 

(37.2) 

179 

(22.4) 

Chemotherapy 

Yes 

28 

(84.9)2 

360 

(47.2) 

22 

(61.1) 

366 

(48.2) 

No 

5 

(15.2) 

402 

(52.8) 

14 

(38.9) 

393 

(51.8) 

Hormonal 

Yes 

21 

(63.6)2 

629 

(83.1) 

25 

(67.6)1 

625 

(83.0) 

treatment 

No 

12 

(36.4) 

128 

(16.9) 

12 

(32.4) 

128 

(17.0) 

171  participants  were  missing  questionnaire 
Fisher’s  exact  tests,  as  appropriate;  ^Among 


data;  ^P<0.05,  2p<0.01,  3p<0.0001, 
postmenopausal  women  only. 


using  chi-squared  o 


AG/GG/AG+GG  AG/AA/AG+AA  AG/AG+AA  AG  GA/AA/GA+AA  TC/TT/TC+TT  AG/AA/AG+AA 


'Figures  1  and  2:  No  participants  had  the  homozygous  variant  genotype  for  SIPA1  rs75894763.  There 
were  no  events  among  participants  with  the  homozygous  variant  genotype  of  SIPA1  r83741376  and 
BRMS1  rs1052S66  in  the  TTR  and  OS  analyses,  respectively. 

Participants  with  the  homozygous  wildtype  genotype  were  used  as  the  reference  group  (not  shown). 
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Abstract  Lymph  node  metastases  and  tumor  character¬ 
istics  predict  breast  cancer  prognosis  but  correlate  imper¬ 
fectly  with  likelihood  of  metastatic  relapse.  Discovery  of 
genetic  polymorphisms  affecting  metastasis  may  improve 
identification  of  patients  requiring  aggressive  adjuvant 
therapy  to  prevent  recurrence.  We  investigated  associations 
between  several  variants  in  the  BRMSl  and  SIPAl  metas¬ 
tasis-modifying  genes  and  lymph  node  metastases,  tumor 
subtype  and  grade,  recurrence,  disease-free  survival,  and 
overall  survival.  This  cross-sectional  and  prospective 
prognostic  analysis  included  859  patients  who  received 
surgery  for  incident  breast  cancer  at  Roswell  Park  Cancer 
Institute,  participated  in  the  DataBank  and  BioRepository 
shared  resource,  and  had  DNA,  clinical,  and  pathology 
data  available  for  analysis.  Genotyping  for  BRMSl 
(rsl  1537993,  rs31 16068,  and  rsl052566)  and  SIPAl 
(rs75894763,  rs746429,  rs3741378,  and  rs2306364)  poly¬ 
morphisms  was  performed  using  Sequenom®  iPLEX  Gold 
and  Taqman®  real-time  PCR  assays.  Logistic  and  Cox 
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proportional  hazards  regressions  were  used  to  estimate 
odds  ratios  (OR)  and  hazard  ratios  (HR),  respectively. 
BRMSl  rsl052566  heterozygous  individuals  were  more 
likely  to  have  node-positive  tumors  (OR  =  1.58,  95  %  Cl 
1.13-2.23),  although  there  was  no  dose-response  rela¬ 
tionship,  and  those  with  at  least  one  variant  allele  were  less 
likely  to  have  the  luminal  B  subtype  (AG  -I-  AA: 
OR  =  0.59,  95  %  Cl  0.36-0.98).  BRMSl  rs31 16068  was 
associated  with  increased  likelihood  of  having  the  luminal 
B  and  the  HER2-enriched  tumor  subtype  (Ptrend  =  0.03). 
Two  SIPAl  SNPs,  rs746429  and  rs2306364,  were  associ¬ 
ated  with  decreased  risk  of  triple-negative  tumors 
(T’trend  =  0.04  and  0.07,  respectively).  Presence  of  8  or 
more  risk  alleles  was  associated  with  an  increased  likeli¬ 
hood  of  having  a  node-positive  tumor  (OR  =  2.14,  95  % 
Cl  1.18-3.36,  Ptrend  =  0.002).  There  were  no  significant 
associations  with  survival.  Polymorphisms  in  metastasis- 
associated  genes  may  be  related  to  tumor  characteristics 
and  lymph  node  metastasis,  but  not  survival.  Future 
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evaluation  of  metastasis-modifying  gene  variants  is  nec¬ 
essary  to  better  understand  the  biology  of  metastasis. 

Keywords  Breast  cancer  •  Metastasis  •  Single  nucleotide 
polymorphism  •  Recurrence  •  Survival 

Abbreviations 

BMI  Body  mass  index 

BRMSl  Breast  cancer  metastasis  suppressor  1 

Cl  Conhdence  interval 

DBBR  DataBank  and  BioRepository 

DCIS  Ductal  carcinoma  in  situ 

DFS  Disease-free  survival 

ER  Estrogen  receptor 

HER2  Human  epidermal  growth  factor  receptor  2 

HR  Hazard  ratio 

HRT  Hormone  replacement  therapy 

NCCN  National  Comprehensive  Cancer  Network 

OR  Odds  ratio 

OS  Overall  survival 

PR  Progesterone  receptor 

RPCl  Roswell  Park  Cancer  Institute 

SlPAl  Signal-induced  proliferation-associated  1 

SNP  Single  nucleotide  polymorphism 

TTR  Time  to  recurrence 


Introduction 

While  early  stage  breast  cancer  has  excellent  prognosis,  it 
is  incurable  once  distant  metastasis  has  occurred  [1,  2]. 
Metastasis  to  regional  lymph  nodes  is  correlated  with  a 
higher  risk  of  developing  distant  metastases  [3,  4],  as  are 
tumor  size  and  grade,  estrogen  and  progesterone  receptor 
status  (ER  and  PR,  respectively),  and  HER2  amplihcation. 
In  general,  larger  tumors  are  correlated  with  increased 
likelihood  of  lymph  node  metastases  at  diagnosis  and  dis¬ 
tant  metastases  [5,  6],  while  ER,  PR,  and  HER2  status  are 
markers  of  tumor  aggressiveness  and  also  determine  suit¬ 
ability  for  targeted  treatments  [7]. 

Even  with  these  known  prognostic  factors,  however,  the 
patients  who  will  ultimately  experience  a  recurrence  are 
not  clearly  identified.  Genetic  variability  may  explain  some 
of  this  heterogeneity  in  metastatic  ability,  particularly  in 
genes  affecting  the  metastatic  cascade.  Many  metastasis- 
related  genes  have  been  identified,  two  of  which  are 
BRMSl  (breast  cancer  metastasis  suppressor  1)  and  SIPAl 
(signal-induced  proliferation-associated  1). 

BRMSl  can  function  as  a  metastasis  suppressor  gene 
[8-10]  that  affects  apoptosis,  colonization,  cell  adhesion, 
and  invasive  potential  by  mitigating  the  effects  of  anti- 
apoptotic  gene  products  regulated  by  the  NF-kB  pathway 


[11,  12].  No  studies  examining  single  nucleotide  poly¬ 
morphisms  (SNPs)  in  BRMSl  have  been  published, 
although  several  expression  studies  have  analyzed  the 
relationship  between  BRMSl  and  breast  tumor  character¬ 
istics  and  prognosis  [13-17].  SIPAl  can  affect  metastatic 
efficiency  by  modifying  cell  adhesion  [18]  and  expression 
of  extracellular  matrix  genes  [19]  and  has  been  shown  to 
promote  metastasis  in  vivo  [20].  Several  SIPAl  SNP 
association  studies  have  been  published,  with  conflicting 
reports  with  respect  to  prognosis  [21-23].  These  data 
indicate  that  BRMSl  and  SIPAl  abnormalities  could  affect 
tumor  aggressiveness,  metastasis,  and  the  risk  of  recurrence 
in  breast  cancer  patients. 

Based  on  this  previously  published  data,  we  selected 
several  SNPs  in  BRMSl  and  SIPAl  to  investigate  as 
potential  candidate  markers  of  tumor  aggressiveness  and 
recurrence.  To  investigate  these  relationships,  we  analyzed 
three  SNPs  mBRMSl  [rsll537993  (Leu67Leu);  rs3116068 
(3'  UTR);  and  rsl052566  (Ala273Val)]  and  four  SNPs  in 
SIPAl  [rs75894763  (Val621Val);  rs746429  (Ala920Ala); 
rs3741378  (Serl82Phe);  and  rs2306364  (Ala342Ala)]  with 
respect  to  lymph  node  metastasis,  tumor  grade  and  subtype, 
time  to  recurrence,  disease-free  survival,  and  overall  sur¬ 
vival  in  women  diagnosed  with  primary,  incident  breast 
cancer. 


Methods 

Study  population  and  outcomes 

We  identified  859  women  diagnosed  between  October 
2003  and  May  2010  with  stage  I-III  incident,  primary, 
histologically  confirmed  breast  cancer,  who  received  sur¬ 
gery  and  treatment  at  Roswell  Park  Cancer  Institute 
(RPCI),  provided  informed  consent  to  RPCFs  DataBank 
and  BioRepository  (DBBR),  and  had  a  DNA  sample 
available.  The  DBBR,  as  previously  described  [24],  is  a 
comprehensive  data  and  sample  bank  containing  high- 
quality  pre-treatment  biospecimens  and  associated  clinical 
and  epidemiologic  data.  All  patients  diagnosed  with  cancer 
at  RPCI  are  invited  to  participate.  After  consent  and  prior 
to  treatment,  including  surgery,  blood  samples  are  col¬ 
lected,  processed,  and  aliquoted  for  storage  in  liquid 
nitrogen.  Epidemiologic  data  obtained  via  self-adminis¬ 
tered  questionnaire  were  available  for  688  of  the  partici¬ 
pants  in  this  analysis. 

Outcomes  were  lymph  node  metastases,  tumor  subtype, 
tumor  grade,  time  to  recurrence,  disease-free  survival,  and 
overall  survival.  Time  to  recurrence  was  defined  as  the 
time  from  diagnosis  to  date  of  first  recurrence  (local  and 
regional  recurrence  and  development  of  distant  metastases) 
or  last  follow-up.  Disease-free  survival  was  defined  as  the 
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time  from  diagnosis  to  the  date  of  first  recurrence,  death 
from  any  cause,  or  last  follow-up.  Overall  survival  was 
defined  as  the  time  from  diagnosis  to  the  date  of  death  from 
any  cause  or  last  follow-up.  Clinical  data  were  obtained 
from  RPCI  clinical  databases  and  supplemented  with  data 
abstracted  from  medical  records  and  the  RPCI  Tumor 
Registry.  Vital  status  and  recurrence  data  were  obtained 
from  the  RPCI  Tumor  Registry  and  the  National  Com¬ 
prehensive  Cancer  Network  (NCCN)  Breast  Cancer  Out¬ 
comes  Database.  The  RPCI  Tumor  Registry  conducts 
yearly  follow-up  on  patients  who  were  last  seen  at  RPCI 
13  months  prior  and  known  to  be  alive.  Vital  status  and 
recurrences  are  ascertained  via  RPCI  medical  record 
abstraction.  Social  Security  Death  Index  and  Legacy.com 
searches,  and/or  letters  sent  to  the  patient,  the  patient’s 
physician,  or  a  family  member.  NCCN-coordinated  linkage 
with  the  National  Death  Index  for  patients  defined  as  “lost 
to  follow-up”  was  completed  on  8  December  2011. 

Fifteen  participants  missing  HER2  status  and  1  missing 
ER  and  PR  status  could  not  be  classified  by  subtype. 
Additionally,  10  participants  were  missing  lymph  node 
status  and  15  were  missing  tumor  grade.  Vital  status  was 
available  for  all  participants.  Eollow-up  ended  in  July  2012 
and  follow-up  time  ranged  from  4  to  101  months.  This 
study  was  approved  by  the  RPCI  Institutional  Review 
Board. 

Genotyping 

The  NCBI  dbSNP  resource  was  used  to  identify  SNPs  in 
the  BRMSl  and  SIPAl  genes  [25].  We  initially  selected  13 
SNPs  In  BRMSl  [rsl7850564  (Aspl75Asp);  rsl  1537993 
(Leu67Leu);  rs75053504  (A/G);  rs31 16068  (A/G);  and 
rsl052566  (C/T))  and  SIPAl  (rs3741378  (Serl82Phe); 
rs76570058  (Prol038Thr);  rs75861149  (Gly368Gly); 
rs2306364  (Ala342Ala);  rs75894763  (Val621Val); 
rs746429  (Ala920Ala);  rs77600626  (Gly249Glu);  and 
rs76089059  (Ala997Ala)]  for  genotyping,  based  on  pres¬ 
ence  in  protein  coding,  3'  untranslated,  or  promoter 
regions,  and  heterozygosity  of  >0.10. 

Genotyping  of  all  13  SNPs  was  conducted  by  the  RPCI 
Genomics  Eacility  using  Sequenom  MassARRAY®  iPLEX 
Gold  matrix-assisted  laser  desorption-ionization  time-of- 
flight  mass  spectrometry  assays.  Genotyping  of  several 
SNPs  (rsl052566,  rs746429,  rs3471378,  rs2306364, 
rs77600626,  and  rs76089059)  was  unsuccessful  using  this 
platform,  and  an  additional  four  SNPs  (rsl7850564, 
rs75053504,  rs75861149,  and  rs76570058)  were  mono- 
morphic  and  not  analyzed  further.  Probes  for  Taqman® 
(Applied  Biosciences)  real-time  PCR  genotyping  assays 
were  available  for  four  of  the  SNPs  that  failed  Sequenom® 
genotyping  (rsl052566,  rs746429,  rs3471378,  and 

rs2306364).  Therefore,  we  were  ultimately  able  to  obtain 


genotyping  data  for  analysis  of  seven  SNPs  using  either  the 
Sequenom®  (rsll537993,  rs3116068,  rs75894763)  or 
Taqman®  (rsl052566,  rs746429,  rs3471378,  and 

rs2306364)  platforms.  Two  SNPs  in  SIPAl,  rs746429,  and 
rs2306364  were  in  strong  linkage  disequilibrium 
(r^  =  0.809).  Duplicate  samples  were  genotyped  to  assess 
intra-  and  inter-plate  reliability.  Genotyping  call  rates 
ranged  from  96.2  to  99.8  %. 

Cross-sectional  analysis 

Statistical  analyses  were  performed  using  SAS  version  9.3. 
Demographic  variables  and  tumor  characteristics  were 
compared  by  lymph  node  status,  tumor  subtype,  and  tumor 
grade  using  Chi-squared  and  Eisher’s  exact  tests  as 
appropriate.  Complete-case  regression  techniques  were 
used  to  analyze  the  relationships  between  SNPs  and  lymph 
node  metastasis,  tumor  grade,  and  tumor  subtype.  Potential 
covarlates  Included  age  at  diagnosis,  tumor  size,  tumor 
grade,  ER,  PR,  and  HER2  status,  lymph  node  metastasis, 
race,  education,  body  mass  index  (BMI),  age  at  menarche, 
menopausal  status,  age  at  menopause,  parity,  age  at  first 
birth,  family  history  of  breast  cancer,  history  of  benign 
breast  disease,  and  hormone  replacement  therapy  (HRT) 
use.  Primary  analyses  incorporated  data  from  the  entire 
study  population  of  859  participants.  Because  we  con¬ 
ducted  complete-case  analyses,  participants  missing  epi¬ 
demiologic  questionnaire  data  dropped  out  of  models 
including  epidemiologic  covariates.  To  minimize  bias 
potentially  introduced  by  these  missing  data,  we  initially 
included  only  tumor  characteristics  in  adjusted  models,  as 
these  data  were  available  for  the  majority  of  our  sample. 
Epidemiologic  variables  were  included  in  separate  models 
to  assess  the  effect  of  their  inclusion  on  odds  ratio  esti¬ 
mates.  A  variable  was  included  in  adjusted  models  if  it  was 
associated  with  the  outcome  and/or  SNP(s),  using  Chi- 
squared  or  Fisher’s  exact  tests  of  significance.  Participants 
with  and  without  questionnaire  data  had  similar  distribu¬ 
tions  of  tumor  characteristics  and  treatment  modalities. 
Participants  missing  questionnaires  were  slightly  younger 
and  somewhat  more  likely  to  have  node-positive  tumors, 
but  differences  were  not  significant.  For  all  outcomes, 
sensitivity  analyses  in  which  we  excluded  participants  who 
self-identified  as  non- white  (5.4  %)  were  performed. 

Unconditional  logistic  regression  was  used  to  estimate 
odds  ratios  (OR)  and  95  %  confidence  intervals  (Cl)  for  the 
associations  between  each  of  the  seven  SNPs  and  lymph 
node  status  and  tumor  grade.  'We  first  constructed  age- 
adjusted  models  and  then  subsequently  added  ER  and  PR 
status,  tumor  size,  and  tumor  grade.  A  third  model  addi¬ 
tionally  included  HER2  status,  race,  education,  HRT,  and 
menopausal  status.  Einally,  we  restricted  our  analysis  to 
include  only  stage  2  and  3  participants,  as  these  patients  are 
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eligible  to  have  lymph  node  metastases  (by  definition,  stage 
1  is  node  negative). 

Moderate-  and  low-grade  tumors  were  combined,  cre¬ 
ating  a  dichotomous  grade  variable  with  categories  of  low/ 
moderate  grade  (well-differentiated  and  moderately  dif¬ 
ferentiated  tumors)  and  high  grade  (poorly  differentiated 
and  undifferentiated  tumors).  Using  a  similar  strategy  as 
outlined  above,  we  first  adjusted  for  age  and  tumor  size.  In 
separate  models,  we  added  nodal  status,  ER,  PR,  and  HER2 
status,  race,  HRT,  and  menopausal  status. 

Generalized  logit  multinomial  logistic  regression  was 
used  to  examine  associations  between  each  SNP  and  tumor 
subtype,  using  the  luminal  A  subtype  as  the  comparison 
group.  Adjusted  models  included  age  and  tumor  size;  age, 
tumor  size,  and  nodal  status;  and  age,  tumor  size,  nodal 
status,  HRT,  race,  and  menopausal  status. 

Eor  all  analyses,  Ptrend  was  calculated  by  coding  geno¬ 
types  as  0,  1,  or  2  for  homozygous  wild-type,  heterozygous, 
or  homozygous  variant  genotypes,  respectively,  and  treating 
the  SNP  as  a  continuous  variable  in  the  regression  models. 

Survival  analysis 

Cox  proportional  hazards  regression  was  used  to  estimate 
hazard  ratios  (HR)  and  95  %  confidence  intervals  (Cl)  for 
the  relationships  between  the  SNPs  and  overall  survival, 
time  to  recurrence,  and  disease-free  survival.  Log-rank 
tests  were  used  to  identify  predictors  for  inclusion  in 
multivariate  proportional  hazards  regression  models. 
Variables  tested  as  predictors  were  age  at  diagnosis,  tumor 
size,  tumor  grade,  ER,  PR,  and  HER2  status,  race,  educa¬ 
tion,  BMI,  age  at  menarche,  menopausal  status,  parity, 
family  history  of  breast  cancer,  history  of  benign  breast 
disease,  hormone  replacement  therapy  use,  radiation 
treatment,  chemotherapy,  hormonal  treatment  (tamoxifen, 
etc.),  and  Charlson  Comorbidity  Score. 

Significant  variables  were  age,  ER  and  PR  status,  tumor 
grade,  tumor  size,  comorbidity  score,  radiation  treatment, 
chemotherapy,  hormone  treatment,  education,  age  at 
menarche,  and  menopausal  status.  Between  7  and  12 
events,  depending  on  the  outcome,  occurred  among  par¬ 
ticipants  missing  chemotherapy,  hormone  treatment,  or 
comorbidity  score.  Adjustment  for  these  variables  would 
therefore  result  in  loss  of  a  large  number  of  events.  Simi¬ 
larly,  adjustment  for  epidemiologic  variables  resulted  in  a 
large  proportion  of  participants  dropping  out  of  our  anal¬ 
yses  due  to  missing  data.  To  minimize  bias  due  to  dropout 
and  maximize  power,  we  initially  limited  model  covariates 
to  age,  ER,  PR,  tumor  grade,  tumor  size,  and  radiation 
treatment.  We  then  tested  the  addition  of  other  covariates 
(comorbidity  score,  chemotherapy,  hormone  therapy,  edu¬ 
cation,  age  at  menarche,  and  menopausal  status)  to  assess 
their  impact  on  the  SNP-survival  outcome  associations. 


Finally,  we  conducted  sensitivity  analyses  in  models 
adjusted  for  age,  ER,  PR,  radiation,  tumor  grade,  and  tumor 
size  by  testing  the  effect  of  including  BMI  and  excluding 
non-white  participants  on  hazard  ratio  estimates. 

Risk  allele  score  construction 

We  constructed  a  summary  risk  allele  score  using  the  log- 
additive  model  to  estimate  unadjusted  per  copy  variant 
allele  odds  ratios  and  hazard  ratios  for  lymph  node  status, 
tumor  grade,  time  to  recurrence,  disease-free  survival,  and 
overall  survival.  We  did  not  include  tumor  subtype  in  this 
analysis  due  to  the  complexity  of  creating  a  summary  score 
for  each  subtype.  Because  the  per  copy  variant  allele 
effects  were  small  in  several  instances,  we  considered  odds 
ratios/hazard  ratios  that  fell  within  the  range  0.95-1.05  as 
being  too  close  to  null  to  assign  a  risk  allele  score.  When 
this  occurred,  0  risk  alleles  were  assigned  for  all  genotypes. 
If  the  per  copy  variant  allele  odds  ratio/hazard  ratio  for  a 
given  SNP  was  greater  than  or  equal  to  1.06,  genotypes  for 
that  SNP  were  assigned  a  score  based  on  the  following 
scheme:  homozygous  wild-type  genotype  =  0  risk  alleles; 
heterozygous  genotype  =  1  risk  allele;  homozygous  vari¬ 
ant  genotype  =  2  risk  alleles.  If  the  odds  ratio/hazard  ratio 
was  less  than  or  equal  to  0.94,  the  coding  scheme  was 
reversed:  homozygous  wild  type  =  2  risk  alleles;  hetero¬ 
zygous  =  1  risk  allele;  homozygous  variant  =  0  risk 
alleles.  The  number  of  risk  alleles  for  each  of  the  seven 
SNPs  was  then  added  together  for  each  participant  to  create 
the  summary  risk  allele  score.  The  directions  of  the  log- 
additive  odds  and  hazard  ratios  for  each  SNP  were  not 
generally  similar  across  the  lymph  node  status,  tumor  grade 
and  survival  analyses,  leading  to  the  assignment  of  dif¬ 
ferent  numbers  of  risk  alleles  for  each  outcome.  Using  the 
distributions  of  risk  alleles  for  the  lymph  node  status,  tumor 
grade,  and  survival  analyses,  we  categorized  the  number  of 
risk  alleles  as  5  or  less,  6,  7,  and  8  or  more  to  create  a 
summary  risk  allele  score.  The  category  “5  or  less”  risk 
alleles  served  as  the  reference  category.  We  then  estimated 
the  odds/hazard  ratios  for  each  level  of  the  risk  allele  score, 
using  logistic  and  proportional  hazards  regression  as 
described  above.  Regression  models  included  the  same 
covariates  as  described  previously  for  the  lymph  node 
status,  tumor  grade,  and  survival  analyses. 

Results 

Participant  characteristics  are  shown  in  Table  1.  Younger 
participants  were  more  likely  to  have  node-positive,  high- 
grade  tumors  of  the  luminal  B,  HER2-enriched,  and  triple¬ 
negative  subtypes.  Higher  educational  attainment  was 
associated  with  decreased  likelihood  of  node-positive 
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tumors.  Ever  users  of  HRT  were  more  likely  to  have 
luminal  B  and  HER2-enriched  tumors.  There  were  no  other 
significant  differences  in  demographic  and  reproductive 
variables  with  respect  to  tumor  characteristics.  Lymph 
node  metastases  were  more  commonly  observed  in  con¬ 
junction  with  high-grade,  larger  size,  ER-/PR-negative,  and 
HER2-positive  tumors. 

SNP  and  lymph  node  status  associations  are  presented  in 
Table  2.  BRMSl  rsl  1537993  and  SIPAl  rs75894763, 
rs746429,  and  rs2306364  were  not  significantly  associated 
with  lymph  node  metastases.  Age-adjusted  odds  ratios 
were  similar  to  those  obtained  in  multivariate  models  (data 
not  shown).  BRMSl  rsl052566  heterozygotes  were  more 
likely  to  have  node-positive  tumors  (OR  =  1.58,  95  %  Cl 
1.13-2.23),  which  remained  significant  after  additional 
adjustment  for  race,  HRT  use,  education,  and  menopausal 
status  (OR  =  1.70,  95  %  Cl  1.13-2.55)  and  when  we 
limited  the  analysis  to  participants  with  stage  2  and  3 
tumors  (OR  =  1.85,  95  %  Cl  1.08-3.18).  However,  this 
relationship  was  not  observed  among  homozygous  indi¬ 
viduals.  Although  only  marginally  significant,  participants 
with  at  least  one  copy  of  the  variant  SIPAl  rs3741378 
allele  were  less  likely  to  have  node-positive  tumors 
(OR  =  0.70,  95  %  Cl  0.48-1.02).  While  the  direction  and 
magnitude  persisted  following  adjustment  for  additional 


covariates  and  limitation  to  stage  2  and  3  tumors,  this 
association  became  nonsignificant.  Similarly,  BRMSl 
rs3116068  approached  statistical  significance  only  when 
race,  HRT  use,  education,  and  menopausal  status  were 
added  as  covariates  (OR  =  0.67,  95  %  Cl  0.45-1.01). 

SNP  associations  with  tumor  grade  are  presented  in 
Table  3.  SIPAl  rs75 894763  heterozygous  participants  were 
more  likely  to  have  high-grade  tumors  (OR  =  2.62,  95  % 
Cl  1.06-6.48),  but  only  after  further  adjustment  for  race, 
HRT  use,  and  menopausal  status. 

Associations  between  tumor  subtype  and  genotype  are 
shown  in  Table  4.  Results  of  age-  and  tumor  size-adjusted 
analyses  were  similar  to  the  findings  presenfed  in  Table  4 
and  are  nol  shown.  The  BRMSl  rs3 116068  homozygous 
varianf  genotype  was  associated  with  increased  likelihood 
of  luminal  B  tumors  (OR  =  2.50,  95  %  Cl  1.10-5.66), 
although  there  was  a  nonsignificant  inverse  relationship 
among  heterozygotes.  Those  who  were  heterozygous  or 
homozygous  were  also  more  likely  to  have  tumors  of  the 
HER2-enriched  subtype  (OR  =  2.45,  95  %  Cl  1.18-5.06, 
T’trend  =  0.03).  These  relationships  remained  significant 
after  additional  adjustment  for  race,  HRT  use,  and  meno¬ 
pausal  status  (data  not  shown).  Patients  homozygous  or 
heterozygous  for  BRMSl  rs  1052566  were  less  likely  to 
have  luminal  B  tumors  (OR  =  0.59,  95  %  Cl  0.36-0.98, 


Table  1  Participant  characteristics  by  lymph  node  status,  tumor  subtype,  and  tumor  grade 


Characteristic^, 
n  (%) 

Lymph  node  status^ 

Tumor  subtype'’ 

Tumor  grade 

b 

Positive 
(N  =  246) 

Negative 
{N  =  603) 

Luminal  A 
{N  =  596) 

Luminal  B 
(N  =  75) 

HER2  (-I-) 
(N  =  34) 

Triple  (— ) 
{N=  138) 

Low 

{N  =  225) 

Moderate 
{N=  371) 

High 

{N  =  248) 

Age  at  diagnosis 

<50 

97  (39.4) 

111  (29.4)* 

174  (29.2) 

29  (38.7) 

15  (44.1) 

51  (37.0) 

62  (27.6) 

115  (31.0) 

93  (37.5) 

51-65 

89  (36.2) 

253  (42.0) 

246  (41.3) 

28  (37.3) 

15  (44.1) 

50  (36.2) 

95  (42.2) 

141  (38.0) 

99  (39.9) 

>66 

60  (24.4) 

173  (28.7) 

176  (29.5) 

18  (24.0) 

4  (11.8) 

37  (26.8) 

68  (30.2) 

115  (31.0) 

56  (22.6) 

Race 

White 

171  (90.5) 

463  (94.5) 

450  (94.5) 

60  (92.3) 

19  (86.4) 

103  (91.2) 

172  (96.1) 

277  (93.9) 

184  (90.6) 

Non-white 

18  (9.5) 

27  (5.5) 

26  (5.5) 

5  (7.7) 

3  (13.6) 

10  (8.8) 

7  (3.9) 

18  (6.1) 

19  (9.4) 

Missing 

57 

113 

120 

10 

12 

25 

46 

76 

45 

Education 

High  school  or  less 

79  (41.8) 

166  (34.5) 

180  (38.4) 

19  (29.7) 

8  (36.4) 

39  (34.8) 

63  (36.0) 

109  (37.3) 

76  (37.8) 

At  least  some  college 

110  (58.2) 

315  (65.5) 

289  (61.6) 

45  (70.3) 

14  (63.6) 

73  (65.2) 

112  (64.0) 

183  (62.7) 

125  (62.2) 

Missing 

57 

122 

127 

11 

12 

26 

50 

79 

47 

Menopausal  status 

Pre 

72  (38.7) 

157  (32.2) 

150  (31.8) 

29  (44.6) 

7  (31.8) 

39  (34.8) 

58  (32.8) 

92  (31.4) 

77  (38.1) 

Post 

114  (61.3) 

331  (67.8) 

322  (68.2) 

36  (55.4) 

15  (68.2) 

73  (65.2) 

119  (67.2) 

201  (68.6) 

125  (61.9) 

Missing 

60 

115 

124 

10 

12 

26 

48 

78 

46 

HRT  use' 

Never 

68  (60.7) 

170  (53.1) 

183  (58.8) 

15  (42.9) 

4  (26.7) 

39  (54.2)* 

66  (57.4) 

113  (58.2) 

62  (50.4) 

Ever 

44  (39.3) 

150  (46.9) 

128  (41.2) 

20  (57.1) 

11  (73.3) 

33  (45.8) 

49  (42.6) 

81  (41.8) 

61  (49.6) 

Missing 

2 

11 

11 

1 

0 

1 

4 

7 

2 

Parity 

Nulliparous 

28  (15.1) 

94  (19.2) 

84  (17.8) 

17  (26.2) 

3  (13.6) 

14  (12.5) 

34  (19.0) 

51  (17.4) 

35  (17.4) 
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Table  1  continued 

Characteristic'*, 

n  (%) 

Lymph  node  status^ 

Tumor  subtype'^ 

Tumor  grade 

b 

Positive 
{N  =  246) 

Negative 
{N  =  603) 

Luminal  A 
{N  =  596) 

Luminal  B 

(N  =  75) 

HER2  (+) 
(N  =  34) 

Triple  (-) 
(N  =  138) 

Low 

{N  =  225) 

Moderate 
{N=  371) 

High 

(N  =  248) 

Parous 

158  (84.9) 

395  (80.8) 

389  (82.2) 

48  (73.8) 

19  (86.4) 

98  (87.5) 

145  (81.0) 

242  (82.6) 

166  (82.6) 

Missing 

60 

114 

123 

10 

12 

26 

46 

78 

47 

ER  status 

Positive 

185  (75.2) 

479  (80.0) 

590  (99.0) 

73  (97.3) 

0 

0 

220  (98.2) 

332  (89.5) 

113  (45.7)* 

Negative 

61  (24.8) 

120  (20.0) 

6  (1.0) 

2  (2.7) 

34 

138 

4  (1.8) 

39  (10.5) 

134  (54.3) 

Missing 

0 

4 

0 

0 

0 

0 

1 

0 

1 

PR  status 

Positive 

154  (62.6) 

420  (70.1)* 

512  (85.9) 

63  (84.0) 

0 

0 

200  (89.3) 

281  (75.7) 

91  (36.8)* 

Negative 

92  (37.4) 

179  (29.9) 

84  (14.1) 

12  (16.0) 

34 

138 

24  (10.7) 

90  (24.3) 

156  (63.2) 

Missing 

0 

4 

0 

0 

0 

0 

1 

0 

1 

HER2  status 

Positive 

44  (17.9) 

65  (11.0)* 

0 

75 

34 

0 

4  (1.8) 

46  (12.6) 

55  (22.4)* 

Negative 

202  (82.1) 

524  (89.0) 

596 

0 

0 

138 

217  (98.2) 

320  (87.4) 

190  (77.6) 

Missing 

0 

14  (2.3) 

0 

0 

0 

0 

4 

5 

3 

Lymph  node  status 

Positive 

246 

- 

163  (27.7) 

23  (30.7) 

21  (61.8) 

39  (28.5)* 

47  (21.4) 

106  (28.8) 

90  (36.6)* 

Negative 

- 

603 

425  (72.3) 

52  (69.3) 

13  (38.2) 

98  (71.5) 

173  (78.6) 

262  (71.2) 

156  (63.4) 

Missing 

- 

- 

8 

0 

0 

1 

5 

3 

2 

Tumor  subtype 

Luminal  A 

163  (66.3) 

425  (72.3)* 

596 

- 

- 

- 

213  (96.8) 

293  (80.1) 

84  (34.3)* 

Luminal  B 

23  (9.4) 

52  (8.8) 

- 

75 

- 

- 

4  (1.8) 

37  (10.1) 

32  (13.1) 

HER2  (+) 

21  (8.5) 

13  (2.2) 

- 

- 

34 

- 

0 

9  (2.5) 

23  (9.4) 

Triple  negative 

39  (15.9) 

98  (16.7) 

- 

- 

- 

138 

3  (1.4) 

27  (7.4) 

106  (43.3) 

Missing 

0 

15 

- 

- 

- 

- 

5 

5 

3 

Tumor  grade 

Low 

47  (19.3) 

173  (29.3)* 

213  (36.1) 

4  (5.5) 

0 

3  (2.2)* 

225 

- 

- 

Moderate 

106  (43.6) 

262  (44.3) 

293  (49.7) 

37  (50.7) 

9  (28.1) 

27  (19.9) 

- 

371 

- 

High 

90  (37.1) 

156  (26.4) 

84  (14.2) 

32  (43.8) 

23  (71.9) 

106  (77.9) 

- 

- 

248 

Missing 

3 

12 

6 

2 

2 

2 

- 

- 

- 

Tumor  size 

Tmi,  TIA  (<5  mm) 

14  (5.7) 

89  (14.8)* 

67  (11.2) 

6  (8.0) 

5  (14.7) 

14  (10.1)* 

37  (16.4) 

40  (10.8) 

18  (7.3)* 

TIB,  TIC  (>5-20  mm) 

116  (47.2) 

413  (68.5) 

398  (66.8) 

45  (60.0) 

14  (41.2) 

78  (56.5) 

159  (70.7) 

239  (64.4) 

136  (54.8) 

T2  (>20-50  mm) 

96  (39.0) 

96  (15.9) 

115  (19.3) 

23  (30.7) 

12  (35.3) 

41  (29.7) 

24  (10.7) 

82  (22.1) 

84  (33.9) 

T3  (>50  mm) 

20  (8.1) 

5  (0.8) 

16  (2.7) 

1  (1.3) 

3  (8.8) 

5  (3.6) 

5  (2.2) 

10  (2.7) 

10  (4.0) 

Stage 

1 

2  (0.9) 

498  (83.4) 

366  (62.8) 

42  (56.0) 

9  (26.5) 

71  (51.8)* 

156  (72.2) 

225  (61.3) 

111  (45.1)* 

2 

164  (66.9) 

99  (16.6) 

170  (29.2) 

27  (36.0) 

15  (44.1) 

50  (36.5) 

48  (22.2) 

115  (31.3) 

96  (39.0) 

3 

79  (32.2) 

0 

47  (8.0) 

6  (8.0) 

10  (29.4) 

16  (11.7) 

12  (5.6) 

27  (7.4) 

39  (15.9) 

Missing 

1 

6 

13 

0 

0 

1 

9 

4 

2 

HRT  hormone  replacement  therapy,  ER  estrogen  receptor,  PR  progesterone  receptor 

*  Significant  at  a  =  0.05;  p  values  were  obtained  from  Chi-squared  and  Fisher’s  exact  tests  as  appropriate  (missing  categories  are  excluded  from p  value 
calculations) 

^  Race,  education,  menopausal  status,  HRT  use,  and  parity  were  available  for  688  participants;  age  at  diagnosis,  ER,  PR,  HER2  status,  tumor  grade,  and 
tumor  size  were  available  for  859  participants 

^  Lymph  node  status,  subtype,  and  grade  were  missing  for  n  =  10,  n  =  16,  and  n  =  15,  respectively 
Excludes  premenopausal  women 
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Table  2  Associations  of  BRMSl  and  SIPAl  SNPs  with  presence  of  lymph  node  metastases  at  diagnosis 


Genotype 

Odds  ratio  for  likelihood  of  node-positive  tumor  at  diagnosis 

Node 
positive,  n 

Node 

negative,  n 

Adjusted”  OR 
(95  %  Cl) 

Adjusted'^  OR 
(95  %  Cl) 

Adjusted”  OR 
(95  %  Cl) 

BRMSl 

rsl  1537993 

AA 

125 

293 

1.00 

1.00 

1.00 

(Leu67Leu) 

AG 

98 

243 

0.93  (0.66-1.30) 

1.02  (0.69-1.53) 

0.98  (0.60-1.63) 

GG 

20 

50 

0.90  (0.49-1.63) 

1.01  (0.53-1.92) 

1.14  (0.44-2.93) 

AG  -f  GG 

0.92  (0.67-1.27) 

1.02  (0.70-1.49) 

1.01  (0.62-1.63) 

^trend 

0.61 

0.95 

0.89 

rs31 16068 

GG 

165 

365 

1.00 

1.00 

1.00 

(3'  UTR) 

AG 

66 

183 

0.80  (0.56-1.14) 

0.66  (0.43-1.02) 

0.69  (0.41-1.18) 

AA 

10 

34 

0.66  (0.30-1.43) 

0.73  (0.31-1.72) 

0.59  (0.18-1.89) 

AG  -1-  AA 

0.78  (0.55-1.09) 

0.67  (0.45-1.01) 

0.68  (0.41-1.13) 

^trend 

0.13 

0.09 

0.13 

rsl052566 

GG 

102 

285 

1.00 

1.00 

1.00 

(Ala273Val) 

AG 

115 

220 

1.58  (1.13-2.23) 

1.70  (1.13-2.55) 

1.85  (1.08-3.18) 

AA 

21 

61 

0.92  (0.51-1.65) 

1.06  (0.54-2.08) 

0.72  (0.33-1.55) 

AG  -1-  AA 

1.43  (1.03-1.99) 

1.56  (1.06-2.30) 

1.48  (0.91-2.42) 

^trend 

0.25 

0.17 

0.67 

SIPAl 

rs75894763 

GG 

237 

558 

1.00 

1.00 

1.00 

(Val621Val) 

AG 

5 

29 

0.39  (0.14-1.08) 

0.51  (0.16-1.60) 

0.30  (0.09-1.03) 

rs746429 

GG 

110 

267 

1.00 

1.00 

1.00 

(Ala920Ala) 

GA 

102 

239 

0.99  (0.70-1.40) 

1.06  (0.71-1.59) 

0.86  (0.52-1.43) 

AA 

29 

66 

1.09  (0.65-1.83) 

1.16  (0.64-2.12) 

1.31  (0.55-3.14) 

GA  -1-  AA 

1.01  (0.73-1.40) 

1.08  (0.74-1.58) 

0.93  (0.57-1.51) 

^trend 

0.82 

0.61 

0.88 

rs3741378 

cc 

187 

413 

1.00 

1.00 

1.00 

(Serl82Phe) 

TC 

49 

148 

0.71  (0.48-1.04) 

0.75  (0.47-1.19) 

0.73  (0.41-1.29) 

TT 

4 

13 

0.62  (0.18-2.10) 

0.61  (0.14-2.67) 

0.42  (0.09-2.05) 

TC  -1-  TT 

0.70  (0.48-1.02) 

0.74  (0.47-1.16) 

0.70  (0.40-1.21) 

^trend 

0.07 

0.18 

0.15 

rs2306364 

GG 

186 

425 

1.00 

1.00 

1.00 

(Ala342Ala) 

AG 

13 

42 

0.71  (0.35-1.41) 

0.70  (0.34-1.48) 

0.89  (0.30-2.63) 

AA 

42 

99 

0.98  (0.64-1.50) 

0.92  (0.56-1.53) 

1.14  (0.58-2.24) 

AG  -1-  AA 

0.90  (0.62-1.32) 

0.85  (0.55-1.32) 

1.07  (0.59-1.95) 

^trend 

0.78 

0.62 

0.75 

Unconditional  logistic  regression  was  used  to  estimate  odds  ratios  (OR)  and  95  %  confidence  intervals  (Cl)  for  risk  of  node-positive  tumors 
“  Adjusted  for  age  at  diagnosis,  ER,  PR,  tumor  size,  and  tumor  grade 
Adjusted  for  age  at  diagnosis,  ER,  PR,  HER2,  tumor  size,  tumor  grade,  race,  HRT,  education,  and  menopausal  status 

Limited  to  women  with  stage  2  and  3  breast  cancer,  adjusted  for  age  at  diagnosis,  ER/PR  status,  and  tumor  grade.  Tumor  size  was  not  included 
as  a  covariate  due  to  sparse  data 


^trend  =  0.05).  These  relationships  remained  marginally 
significant  after  further  adjustment  for  race,  HRT  use,  and 
menopausal  status  (OR  =  0.59,  95  %  Cl  0.34-1.02).  A 
nonsignificant  decrease  in  the  likelihood  of  the  HER2- 
enriched  suhtype  was  observed  among  patients  heterozy¬ 
gous  or  homozygous  for  this  variant,  which  became  sta¬ 
tistically  significant  when  race,  HRT  use,  and  menopausal 
status  were  included  in  the  model  (OR  =  0.32,  95  %  Cl 
0.12-0.90).  Participants  homozygous  for  SIPAl  rs746429 


were  less  likely  to  have  triple-negative  tumors  (OR  =  0.48, 
95  %  Cl  0.24-0.97,  Ptrenu  =  0.04).  Similarly,  participants 
either  homozygous  or  heterozygous  for  SIPAl  rs2306364 
were  less  likely  to  have  triple-negative  tumors,  although 
this  relationship  was  only  borderline  significant 
(OR  =  0.62,  95  %  Cl  0.38-1.01). 

Associations  with  survival  outcomes  are  shown  in 
Table  5.  The  median  follow-up  time  was  45  months  (range 
4-101  months),  during  which  58  recurrences  and  70  deaths 


^  Springer 


880 


Breast  Cancer  Res  Treat  (2013)  139:873-885 


Table  3 

Associations  of  BRMSl  and  SIPAl  SNPs  with  tumor  grade 

Genotype 

Odds  ratio  for  likelihood  of  high-grade  tumor  at  diagnosis 

High  grade,  n 

Low/moderate 
grade,  n 

Adjusted"*  OR 
(95  %  Cl) 

Adjusted'*  OR 
(95  %  Cl) 

BRMSl 

rsl  1537993 

AA 

128 

286 

1.00 

1.00 

(Leu67Leu) 

AG 

92 

243 

0.86  (0.58-1.27) 

0.97  (0.62-1.53) 

GG 

23 

47 

1.23  (0.65-2.35) 

1.36  (0.68-2.70) 

AG  -f  GG 

0.92  (0.64-1.33) 

1.04  (0.68-1.60) 

f^trend 

0.97 

0.55 

rs3 116068 

GG 

156 

369 

1.00 

1.00 

(3'  UTR) 

AG 

70 

175 

0.83  (0.55-1.26) 

0.98  (0.61-1.56) 

AA 

11 

32 

0.82  (0.36-1.89) 

0.67  (0.26-1.75) 

AG  -1-  AA 

0.83  (0.56-1.23) 

0.92  (0.59-1.44) 

f^trend 

0.38 

0.54 

rsl052566 

GG 

113 

267 

1.00 

1.00 

(Ala273Val) 

AG 

92 

240 

0.86  (0.58-1.28) 

0.97  (0.61-1.53) 

AA 

27 

55 

1.39  (0.76-2.57) 

1.63  (0.82-3.25) 

AG  -1-  AA 

0.95  (0.66-1.38) 

1.08  (0.70-1.66) 

^trend 

0.67 

0.33 

SIPAl 

rs75  894763 

GG 

233 

553 

1.00 

1.00 

(Val621Val) 

AG 

10 

23 

1.97  (0.87^.47) 

2.62  (1.06-6.48) 

rs746429 

GG 

114 

259 

1.00 

1.00 

(Ala920Ala) 

GA 

104 

232 

1.20  (0.81-1.78) 

1.28  (0.82-2.00) 

AA 

21 

73 

0.82  (0.44-1.54) 

0.83  (0.41-1.70) 

GA  -1-  AA 

1.11  (0.77-1.61) 

1.17  (0.77-1.79) 

^trend 

0.98 

0.91 

rs3741378 

cc 

166 

429 

1.00 

1.00 

(Serl82Phe) 

TC 

65 

127 

1.44  (0.94-2.20) 

1.46  (0.89-2.39) 

TT 

6 

11 

0.99  (0.27-3.60) 

1.43  (0.36-5.71) 

TC  -1-  TT 

1.40  (0.92-2.11) 

1.46  (0.90-2.35) 

^trend 

0.18 

0.14 

rs2306364 

GG 

187 

419 

1.00 

1.00 

(Ala342Ala) 

AG 

13 

40 

0.84  (0.40-1.77) 

0.72  (0.32-1.61) 

AA 

36 

103 

0.93  (0.57-1.53) 

0.93  (0.53-1.63) 

AG  -1-  AA 

0.90  (0.58-1.40) 

0.86  (0.52-1.40) 

f^trend 

0.71 

0.67 

Unconditional  logistic  regression  was  used  to  estimate  odds  ratios  (OR)  and  95  %  confidence  intervals  (Cl)  for  risk  of  high-grade  tumors 
Low  grade  =  well  differentiated,  Moderate  grade  =  moderately  differentiated,  High  grade  =  poorly  differentiated  and  undifferentiated 
“  Adjusted  for  age  at  diagnosis,  tumor  size,  lymph  node  status,  ER,  PR,  and  HER2 
Adjusted  for  age  at  diagnosis,  tumor  size,  lymph  node  status,  ER,  PR,  HER2,  race,  HRT,  and  menopausal  status 


from  all  causes  occurred.  When  adjusted  for  age  only,  the 
heterozygous  genotype  of  BRMSl  rs3116068  was  associ¬ 
ated  with  poorer  overall  survival  (HR  =  1.65,  95  %  Cl 
1.02-2.68),  hut  this  association  became  nonsignificant 
when  additional  covariates  were  included.  We  did  not 
observe  any  other  significant  associations.  Results  were 
unchanged  when  additional  covariates  (chemotherapy, 
hormone  therapy,  education,  age  at  menarche,  and  meno¬ 
pausal  status)  were  included  (data  not  shown). 


One  of  our  hypotheses  is  that  these  SNPs  would  be 
related  to  lymph  node  status.  Because  nodal  status  is  also 
related  to  survival  outcomes,  it  can  be  hypothesized  that 
this  variable  is  part  of  the  causal  pathway,  and  therefore 
adjustment  for  nodal  status  could  mask  true  SNP-survival 
associations.  To  test  whether  these  SNPs  could  affect  vital 
status  or  recurrence  through  a  pathway  independent  of 
lymph  node  status,  we  included  nodal  status  in  a  model 
containing  age,  ER,  PR,  radiation,  tumor  grade,  and  tumor 
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Table  4  Associations  of  BRMSl  and  SIPAl  SNPs  with  tumor  subtype 

Genotype  Odds  ratio  for  likelihood  of  luminal  B,  HER2  (+),  and  triple  (— )  subtype,  compared  to  luminal  A 


Luminal 

Luminal 

Adjusted'’  OR 

HER2 

Adjusted'’  OR 

Triple 

Adjusted'’  OR 

A,  n 

B,  n 

(95  %  Cl) 

(+),  n 

(95  %  Cl) 

(-),n 

(95  %  Cl) 

BRMSl 

rsl  1537993 

AA 

292 

33 

1.00 

18 

1.00 

77 

1.00 

(Leu67Leu) 

AG 

246 

31 

1.14 

(0.68-1.93) 

15 

1.13 

(0.55-2.34) 

48 

0.74 

(0.50-1.11) 

GG 

49 

10 

1.85 

(0.85^.01) 

1 

0.36 

(0.05-2.82) 

11 

0.74 

(0.50-1.11) 

AG  -f  GG 

1.26 

(0.77-2.06) 

0.99 

(0.49-2.03) 

0.76 

(0.52-1.12) 

P 

^  trend 

0.17 

0.64 

0.26 

rs3 116068 

GG 

379 

50 

1.00 

15 

1.00 

90 

1.00 

(3'  UTR) 

AG 

174 

16 

0.71 

(0.39-1.28) 

16 

2.49 

(1.18-5.27) 

41 

1.00 

(0.66-1.51) 

AA 

28 

9 

2.50 

(1.10-5.66) 

2 

2.11 

(0.45-10.0) 

5 

0.76 

(0.28-2.04) 

AG  -1-  AA 

0.95 

(0.57-1.59) 

2.45 

(1.18-5.06) 

0.97 

(0.65-1.44) 

P 

^  trend 

0.39 

0.03 

0.74 

rsl052566 

GG 

268 

44 

1.00 

18 

1.00 

58 

1.00 

(Ala273Val) 

AG 

240 

24 

0.61 

(0.36-1.03) 

13 

0.70 

(0.33-1.49) 

57 

1.12 

(0.74-1.68) 

AA 

61 

6 

0.55 

(0.22-1.36) 

1 

0.21 

(0.03-1.66) 

15 

1.07 

(0.57-2.03) 

AG  -1-  AA 

0.59 

(0.36-0.98) 

0.60 

(0.29-1.25) 

1.11 

(0.75-1.63) 

Ptrend 

0.05 

0.09 

0.69 

SIPAl 

rs75894763 

GG 

559 

69 

1.00 

34 

1.00 

135 

1.00 

(Val621Val) 

AG 

27 

5 

1.47 

(0.54-3.99) 

0 

NA 

2 

0.29 

(0.07-1.26) 

rs746429 

GG 

255 

35 

1.00 

19 

1.00 

69 

1.00 

(Ala920Ala) 

GA 

245 

29 

0.85 

(0.50-1.44) 

11 

0.56 

(0.26-1.22) 

55 

0.81 

(0.54-1.20) 

AA 

75 

9 

0.86 

(0.39-1.87) 

3 

0.51 

(0.14-1.80) 

10 

0.48 

(0.24-0.97) 

GA  -1-  AA 

0.85 

(0.52-1.39) 

0.55 

(0.27-1.13) 

0.73 

(0.50-1.07) 

Ptrend 

0.57 

0.13 

0.04 

rs3741378 

cc 

429 

53 

1.00 

27 

1.00 

93 

1.00 

(Serl82Phe) 

TC 

138 

17 

0.99 

(0.55-1.78) 

4 

0.46 

(0.16-1.37) 

37 

1.21 

(0.79-1.87) 

TT 

9 

3 

2.91 

(0.75-11.2) 

2 

5.77 

(1.11-29.8) 

3 

1.57 

(0.41-5.98) 

TC  -1-  TT 

1.11 

(0.64-1.92) 

0.68 

(0.27-1.71) 

1.24 

(0.81-1.88) 

Ptrend 

0.45 

0.90 

0.30 

rs2306364 

GG 

421 

58 

1.00 

27 

1.00 

107 

1.00 

(Ala342Ala) 

AG 

42 

6 

1.02 

(0.41-2.51) 

1 

0.38 

(0.05-2.94) 

6 

0.56 

(0.23-1.36) 

AA 

109 

9 

0.59 

(0.28-1.23) 

5 

0.74 

(0.27-2.00) 

18 

0.64 

(0.37-1.11) 

AG  -1-  AA 

0.71 

(0.39-1.29) 

0.64 

(0.26-1.60) 

0.62 

(0.38-1.01) 

Ptrend 

0.18 

0.44 

0.07 

Multinomial  logistic  regression  (generalized  logit  model)  was  used  to  estimate  odds  ratios  (OR)  and  95  %  confidence  intervals  (Cl)  for  risk  of 
luminal  B,  HER2  (+),  or  triple-negative  tumor  subtype,  using  luminal  A  as  the  comparison  group 

Luminal  A  =  ER  and/or  PR  positive,  HER2  negative;  Luminal  B  =  ER  and/or  PR  positive,  HER2  positive;  HER2-enriched  subtype  =  (HER2  (-f)) 
ER  and  PR  negative,  HER2  positive;  Triple  negative  =  (Triple  (— ))  ER,  PR,  and  HER2  negative 

“  Adjusted  for  age  at  diagnosis,  tumor  size,  and  lymph  node  status 


size  as  covariates.  There  was  no  change  in  the  hazard  ratio 
estimates  for  any  of  the  three  outcomes  (data  not  shown). 

Results  of  the  summary  risk  allele  score  analysis  are 
shown  in  Table  6.  Having  eight  or  more  risk  alleles  was 
associated  with  significantly  greater  likelihood  of  having  a 
node-positive  tumor  (OR  =  2.14  95  %  Cl  1.18-3.86). 
There  was  evidence  of  a  dose-response  pattern,  with 


increasing  numbers  of  risk  alleles  associated  with  increased 
likelihood  of  node  positivity  (Ptrend  =  0.002).  There  were 
no  significant  associations  with  tumor  grade  or  the  three 
survival  outcomes. 

In  sensitivity  analyses,  non-white  participants  were 
excluded  for  all  study  outcomes,  which  did  not  alter  our 
findings  (data  not  shown). 
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Table  5  BRMSl  and  SIPAl  SNP  associations  with  time  to  recurrence. 

disease-free  survival,  and  overall  survival 

Genotype 

No.  events, 
TTR 

TTR,  adjusted^ 
HR  (95  %  Cl) 

No.  events,  DFS 

DFS,  adjusted^ 
HR  (95  %  Cl) 

No.  events, 
OS 

OS,  adjusted” 

HR  (95  %  Cl) 

BRMSl 

rsl  1537993 

AA 

27 

1.00 

39 

1.00 

28 

1.00 

(Leu67Leu) 

AG 

21 

1.05  (0.59-1.89) 

35 

1.12  (0.70-1.79) 

28 

1.10 

(0.64-1.89) 

GG 

5 

1.24  (0.47-3.30) 

7 

1.18  (0.52-2.69) 

6 

1.30 

(0.53-3.20) 

AG  -b  GG 

1.08  (0.62-1.89) 

1.13  (0.72-1.77) 

1.13 

(0.67-1.89) 

^trend 

0.69 

0.59 

0.56 

rs3 116068 

GG 

35 

1.00 

48 

1.00 

33 

1.00 

(3'  UTR) 

AG 

15 

0.77  (0.41-1.43) 

29 

1.04  (0.65-1.67) 

26 

1.44 

(0.85-2.45) 

AA 

3 

0.90  (0.26-3.13) 

3 

0.54  (0.16-1.82) 

1 

0.33 

(0.04-2.49) 

AG  -b  AA 

0.78  (0.43-1.42) 

0.97  (0.61-1.55) 

1.31 

(0.77-2.22) 

^trend 

0.51 

0.60 

0.76 

rsl052566 

GG 

26 

1.00 

41 

1.00 

33 

1.00 

(Ala273Val) 

AG 

22 

1.05  (0.59-1.87) 

33 

0.99  (0.62-1.57) 

23 

0.86 

(0.50-1.48) 

AA 

4 

1.01  (0.35-2.96) 

5 

0.97  (0.38-2.48) 

3 

0.81 

(0.24-2.67) 

AG  -b  AA 

1.04  (0.60-1.81) 

0.98  (0.63-1.54) 

0.85 

(0.51-1.44) 

^trend 

0.92 

0.93 

0.55 

SIPAl 

rs75  894763 

GG 

51 

1.00 

78 

1.00 

60 

1.00 

(Val621Val) 

AG 

2 

1.61  (0.38-6.89) 

3 

1.57  (0.48-5.11) 

2 

1.35 

(0.32-5.75) 

rs746429 

GG 

24 

1.00 

37 

1.00 

28 

1.00 

(Ala920Ala) 

GA 

25 

1.20  (0.67-2.14) 

38 

1.15  (0.72-1.83) 

30 

1.24 

(0.73-2.10) 

AA 

5 

0.94  (0.36-2.50) 

7 

0.86  (0.38-1.94) 

4 

0.69 

(0.24-1.99) 

GA  -b  AA 

1.14  (0.66-1.99) 

1.09  (0.70-1.70) 

1.13 

(0.67-1.89) 

^trend 

0.84 

0.98 

0.94 

rs3741378 

cc 

37 

1.00 

58 

1.00 

45 

1.00 

(Serl82Phe) 

TC 

17 

1.28  (0.71-2.30) 

24 

1.09  (0.67-1.77) 

16 

0.86 

(0.48-1.55) 

TT 

0 

NA 

0 

NA 

1 

1.16 

(0.16-8.63) 

TC  -b  TT 

1.20  (0.67-2.16) 

1.03  (0.63-1.67) 

0.87 

(0.49-1.55) 

^trend 

NA 

NA 

0.70 

rs2306364 

GG 

41 

1.00 

61 

1.00 

48 

1.00 

(Ala342Ala) 

AG 

5 

1.16  (0.45-3.01) 

7 

1.20  (0.54-2.71) 

4 

0.92 

(0.32-2.64) 

AA 

8 

0.82  (0.38-1.75) 

13 

0.91  (0.50-1.66) 

9 

0.80 

(0.39-1.65) 

AG  -b  AA 

0.92  (0.49-1.73) 

0.99  (0.59-1.66) 

0.83 

(0.45-1.56) 

^trend 

0.67 

0.84 

0.54 

Patients  whose  recurrence  status  indicated  that  they  were  never  disease  free  were  excluded  from  time  to  recurrence  and  disease-free  survival 
analyses,  but  were  included  in  overall  survival  analyses 

Hazard  ratios  (HR)  and  95  %  confidence  intervals  (Cl)  were  estimated  using  Cox  proportional  hazards  regression 

TTR  =  Time  to  recurrence  (time  from  diagnosis  to  date  of  first  recurrence  or  date  of  last  follow-up),  DFS  =  disease-free  survival  (time  from 
diagnosis  to  date  of  first  recurrence,  death,  or  last  follow-up),  OS  =  overall  survival  (time  from  diagnosis  to  date  of  death  or  date  of  last  follow¬ 
up) 

“  Adjusted  for  age  at  diagnosis,  ER,  PR,  tumor  size,  tumor  grade,  radiation  treatment,  and  Charlson  Comorbidity  Score 


Discussion 

Our  data  suggest  that  2  SNPs  in  the  BRMSl  gene, 
rsl052566  and  rs3 116068,  may  be  associated  with  lymph 
node  status  and  tumor  subtype,  and  that  SNPs  in  the  SIPAl 
gene  may  be  associated  with  tumor  grade  and  subtype.  We 
also  found  that  a  summary  score,  composed  of  the  number 
of  “at  risk”  alleles  for  each  of  the  seven  BRMSl  and  SIPAl 


SNPs  analyzed,  was  signihcantly  associated  with  lymph 
node  status. 

To  our  knowledge,  associations  between  SNPs  in 
BRMSl  and  breast  tumor  characteristics  and  prognosis 
have  not  been  previously  evaluated.  BRMSl  has  multiple 
functions,  including  transcriptional  regulation  via  NF-kB 
signaling  pathways  [26,  27],  chromatin  modification  [26], 
interactions  with  histone  deacetylase  complexes  [27],  and 
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Table  6  Risk  allele  score  associations  with  lymph  node  status,  tumor  grade,  time  to  recurrence,  disease-free  survival,  and  overall  survival 


Number 
of  alleles 

Lymph  node  status 

Tumor  grade 

Time  to 

recurrence 

Disease-free 

survival 

Overall  survival 

Positive 

(n) 

Negative 

(n) 

OR“ 

(95  %  Cf) 

High 

(n) 

Low/ 

Mod 

(n) 

or'" 

(95  %  Cl) 

N* 

HR" 

(95  %  Cl) 

N* 

HR" 

(95  %  Cl) 

N* 

HR" 

(95  %  Cl) 

5  or  less 

18 

73 

1.00 

153 

362 

1.00 

19 

1.00 

15 

1.00 

34 

1.00 

6 

50 

148 

1.31 

(0.69-2.50) 

53 

150 

0.85 

(0.54-1.32) 

21 

1.83 

(0.96-3.48) 

32 

1.84 

(0.98-3.46) 

20 

1.19 

(0.67-2.10) 

7 

44 

117 

1.38 

(0.71-2.67) 

33 

61 

1.19 

(0.67-2.10) 

12 

1.22 

(0.59-2.54) 

15 

0.86 

(0.41-1.82) 

8 

1.74 

(0.79-3.83) 

8  or  more 

131 

251 

2.14 

(1.18-3.86) 

4 

6 

1.28 

(0.28-5.91) 

2 

0.50 

(0.11-2.17) 

20 

1.13 

(0.57-2.24) 

0 

NA 

P  trend 

0.002 

0.78 

0.82 

0.60 

0.48 

N  is  the  number  of  participants  in  each  risk  allele  score  category 
N*  is  the  number  of  events  in  each  risk  allele  score  category 
“  Adjusted  for  age  at  diagnosis,  ER,  PR,  tumor  size,  and  tumor  grade 
Adjusted  for  age  at  diagnosis,  ER,  PR,  HER2,  tumor  size,  and  lymph  node  status 

Adjusted  for  age  at  diagnosis,  ER,  PR,  tumor  size,  tumor  grade,  radiation  treatment,  and  Charlson  Comorbidity  Score 


transcriptional  repression  of  anti-apoptotic  genes  [12]. 
Previous  studies  have  correlated  decreased  BRMSl  gene 
expression  with  breast  tumor  aggressiveness.  Reduced 
mRNA  and  protein  expression  in  breast  tumors  has  been 
associated  with  PR-negative,  HER2-positive  tumors,  as 
well  as  younger  age  at  diagnosis  [13,  14],  but  not  with 
lymph  node  metastases  [14,  15].  One  study  found  that 
BRMSl  mRNA  was  reduced  in  brain  metastases  of  breast 
cancer  patients  [16].  In  survival  analyses,  both  increased 
and  decreased  BRMSl  gene  expression  have  been  corre¬ 
lated  with  reduced  survival  [14,  17],  while  loss  of  BRMSl 
protein  expression  has  been  associated  with  reduced  sur¬ 
vival  only  in  patients  with  ER-negative,  HER2-positive 
tumors  [13]. 

We  observed  several  relationships  between  SlPAl  SNPs 
rs746429,  rs3741378,  and  rs2306364  and  tumor  subtype, 
although  there  are  too  few  participants  with  variant  alleles 
to  draw  strong  conclusions.  SlPAl  encodes  a  GTPase- 
activating  protein  and  affects  expression  of  extracellular 
matrix  genes  [19].  SlPAl  SNPs  rs931127,  rs3741378,  and 
rs746429  have  not  been  shown  to  be  associated  with 
overall  survival  [22],  similar  to  our  findings.  SlPAl 
rs3741378  has  been  correlated  with  increased  risk  of  ER/ 
PR-negative  tumors,  while  rs746429  has  been  correlated 
with  increased  risk  of  node-positive  breast  tumors  [21].  In 
our  study,  rs374I378  was  associated  with  an  increased  risk 
of  HER2-enriched  tumors,  which  are  ER/PR  negative. 

Our  study  used  data  and  samples  collected  under  the 
DBBR’s  standardized  protocol  and  had  relatively  large 
sample  size,  but  our  analyses  were  limited  by  the  small 
number  of  outcomes  and  participants  with  the  variant 
genotypes.  We  also  examined  only  a  few  SNPs  in  each 


gene.  Our  ability  to  adjust  for  socioeconomic  and  repro¬ 
ductive  covariates  was  limited  by  missing  epidemiologic 
questionnaire  data,  although  these  factors  are  unlikely  to  be 
strong  confounders  as  SNPs  are  generally  unlikely  to  be 
associated  with  them.  Distributions  of  tumor  characteristics 
and  treatment  were  similar  between  participants  with  and 
without  questionnaire  data. 

We  compared  the  characteristics  of  participants  with 
complete  data  to  the  characteristics  of  those  with  incom¬ 
plete  data  to  assess  the  possibility  of  bias.  In  general, 
participants  with  incomplete  data  were  more  likely  to  be 
postmenopausal,  non-white,  to  have  smaller  tumors,  and  to 
have  received  hormonal  treatment,  but  were  less  likely  to 
have  a  family  history  of  breast  cancer  and  to  have  received 
radiation  and  chemotherapy.  Age,  ER,  PR,  HER2,  tumor 
grade,  comorbidity  score,  parity,  age  at  menarche,  and 
history  of  benign  breast  disease  were  generally  similar. 
This  indicates  that  the  participants  who  dropped  out  of  our 
analyses  had  less  aggressive  tumors  than  those  who  were 
included,  but  were  similar  with  respect  to  other  possible 
risk  factors.  Genotype  frequencies  were  also  not  signifi¬ 
cantly  different  between  groups,  suggesting  that  a  sub¬ 
stantial  bias  is  unlikely  to  be  present. 

The  median  follow-up  time  in  this  study  was  45  months 
(3.8  years),  which  may  have  been  too  short  to  observe 
associations  with  recurrence  and  survival  if  there  is  a  true 
effect  of  the  SNPs  on  these  outcomes.  We  did  not  correct 
for  multiple  comparisons  and  have  not  performed  a  repli¬ 
cation  study.  The  types  of  breast  cancer  patients  treated  at 
our  institution  and  community-based  facilities  could  be 
different  with  respect  to  tumor  characteristics,  possibly 
leading  to  our  study  population  having  a  greater  proportion 
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of  aggressive  tumor  characteristics  than  would  be  expected 
from  the  source  population,  which  could  affect  the  gener- 
alizability  of  our  findings.  It  is  possible,  although  less 
likely,  that  there  could  be  differences  with  respect  to 
genotype  as  well. 

Conclusions 

In  conclusion,  we  showed  that  SNPs  in  BRMSl  and  SIPAl 
may  be  associated  with  tumor  characteristics  related  to 
prognosis.  Additional  studies  are  needed  to  validate  these 
findings  and  furfher  investigafe  relationships  between 
genetic  variation  in  metastasis-modifying  genes  and  the 
metastatic  phenotype.  Understanding  the  biology  of 
metastasis  and  identifying  biomarkers  of  recurrence  are 
necessary  to  improve  prediction  of  the  subset  of  patients 
who  will  experience  metastatic  relapse,  particularly  since 
treatment  for  breast  cancer  often  results  in  significant 
patient  morbidity.  While  many  women  will  never  progress 
to  metastatic  disease,  identifying  factors  associated  with 
metastatic  recurrence  is  critical  to  achieving  effective, 
efficient  therapy  for  breast  cancer  patients. 
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